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A Study  of  the  Time-Dependent 

1 

V7ind-Driven  Ocean  Circulation 

by 

2 3 

G.  Veronis  and  G,  W,  Morgan 

Abstract*  This  investigation  is  concerned  with  the 
large-scale  wind-driven  motions  of  the  ocean  and  their  responses 
to  a time  variation  in  the  wind#  Starting  from  the  equations 
of  motion  for  an  inhomogeneous  fluid,  a detailed  formulation  of 
the  problem  is  presented,  Including  the  listing  and  discussion 
of  the  assumptions  and  simplifications  necessary  to  reduce  the 
general  mathematical  model  to  one  which  may  be  successfully 
attacked  analytically# 

Since  the  real  ocean  is  baroclinlQ,  the  problem  is 
formulated  to  include  a non-uniform  density  distribution.  Two 
special  cases  are  considered# 

(i)  An  ocean  consisting  of  two  superposed  layers  of  con- 
stant density  is  assumed  and  the  equations  are  integrated  over 
each  layer  to  simplify  the  analysis*  Attempts  at  an  analytical 
solution  for  this  case  were  unsuccessful# 


^ The  results  presented  in  this  paper  were  obtained  in  the 
course  of  research  conducted  under  Contract  N7onr-35801„ 


p 

Research  Assistant,  Graduate  Division  of  Applied  Mathematics, 
Brown  University,  Providence,  R.  I. 

^ Associate  Professor  of  Applied  Mathematics,  Brown  University, 
Providence,  R.  I. 
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(ii)  A more  general  density  distribution  is  then  assumed, 
but  a more  restrictive  assumption  is  made  concerning  the  verti- 
cal variation  of  velocity.  In  particular,  it  is  assumed  that 
there  exists  a (variable)  depth  below  which  the  velocities  are 
negligible.  As  a result  of  this  assumption,  a direct  relation 
is  found  between  the  thermocline  and  the  free  surface.  The 
equations  are  Integrated  from  this  depth  up  to  the  free  surface. 
The  linearized  equations  are  then  subjected  to  an  analytical 
treatment  consisting  of  a perturbation  expansion  in  terms  of  a 
parameter  which  is  proportional  to  the  frequency  of  the  wind 
variation.  The  resulting  equations  are  solved  by  boundary 
layer  technique. 

Results  are  derived  for  the  response  of  the  mass  trans- 
port to  slowly  varying  winds,  and  the  effect  of  the  wind  on  the 
intensified  stream  near  the  western  boundary  is  discussed  in 
detail* 

The  two-layer  steady  problem  is  also  solved  and  the 
steady  position  of  the  thermocline  is  determined. 
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1,  Introduction.  Much  of  the  investigation,  both  theo- 
retical and  observational  in  the  field  of  oceanography  has  cente]> 
ed  around  the  dynamics  of  ocean  currents  - including  the  mass 
transport  of  the  Gulf  Stream  and  the  Kuroshio  Current,  and  the 
general  oceanic  circulation.  Recently  interest  has  developed 
regarding  the  response  of  the  thermocline  (the  region  of  sharp 
vertical  gradient  of  density)  to  a time-varying  wind. 

Since  the  time  of  Ekman’s  first  paper  [ 1]  j a large 
number  of  papers  have  appeared  in  some  o.f  the  geophysical  jour- 
nals dealing  with  various  aspects  of  ocean  currents®  However, 
analytical  investigacions  of  the  problem  of  general  oceanic 
circulation  have  met  with  success  only  in  recent  years.  In  the 
past  decade  various  Interesting  and  meaningful  mathematical 
models  have  been  suggested  by  numerous  investigators#  Sverdrup 
[2]  and  Reid  [3]  proposed  a fairly  simple  model  which  seems  to 
give  very  good  qualitative  results  for  a region  with  only  one 
north-south  boundary.  Stomniel  [4]  considered  two  linearized 
models  with  a simplified  viscous  term.  His  very  important  con- 
tribution to  the  overall  problem  is  based  on  the  difference 
between  the  results  obtained  with  the  two  models.  In  one  case, 
the  Coriolis  term  was  constant  and  the  resulting  streamline 
pattern  is  identical  with  the  one  in  a model  with  no  rotation. 

In  the  second  case,  the  Coriolis  term  varied  linearly  with 
latitude  and  westward  intensification  resulted  - a factor  which 

♦ Numbers  in  square  brackets  refer  to  the  bibliography  at  the 
end  of  the  paper, 
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was  not  present  in  the  previous  case.  Since  Storamel’s  paper 
all  problems  dealing  with  general  circulation  contain  a varying 
Coriolis  parameter,  Munk  [?]  refined  all  the  previous  work  and 
included  the  general  viscous  terms  in  the  equations  of  motion* 

He  solved  the  problem  of  a steady  wind  blowing  over  an  enclosed 
ocean,  taking  account  of  many  of  the  salient  features  which  are 
present  in  the  real  ocean.  Munk*s  work  was  extended  by  Munk 
and  Carrier  [63  to  Include  oceans  of  various  geometrical  shapes, 
viz* , triangular  and  semi-circular.  It  was  further  extended 
by  Munk, Groves,  and  Carrier  [7]  to  include  the  non-linear  terms 
by  means  of  a perturbation  procedure. 

Along  with  the  American  publications,  a number  of  papers 
have  appeared  in  Japan,  Notable  among  the  Japanese  authors  is 
Hidaka,  who  published  a series  of  articles  covering  many  of  the 
Interesting  phenomena  of  oceanographic  problems,  Among  his  con- 
tributions are  a series  of  three  papers  on  drift  currents  in  an 
enclosed  ocean  [12],  [I3] , [l4],  and  a contribution  concerning 
the  neglect  of  the  non-linear  terms  in  the  solution  of  problems 
in  dynamic  oceanography  [1^3. 

Practically  all  of  the  work  done  so  far  in  ocean  current 
problems  has  been  confined  to  motions  which  are  Independent  of 
time.  Each  publication  has  treated  some  aspect  of  the  general 
problem  of  oceanic  circulation,  This  problem  essentially  con- 
sists of  finding  the  dynamic  pattern  which  results  from  a given 
distribution  of  winds  acting  on  the  ocean  surface. 

The  complete  problem  contains  a large  number  of  features, 
such  as  large-scale  oceanic  circulation,  surface  waves,upwelling. 
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etc.  To  find  all  such  motions  one  would  have  to  take  into  ac- 
count the  effects  of  the  wind,  density  and  temperature  distribu- 
tion, the  topography  of  the  ocean  bed  and  possibly  even  such 
features  as  salinity.  Needless  to  say,  a mathematical  analysis 
Including  all  these  features  is  impossible.  It  is  therefore 
necessary  to  decide  what  particular  aspects  of  the  problem  one 
wishes  to  study.  In  this  paper  we  shall  confine  our  attention 
to  large-scale  wind-driven  motions  in  the  oceans  and  their  re-* 
sponses  to  a prescribed  time  variation  in  the  wind*  In  the 
Atlantic  Ocean,  such  large-scale  motions  must  include  the  Gulf 
Stream  and  its  counter-currents,  the  Sargasso  Sea,  etc. 

The  time-dependent  problem  has  also  been  considered  by 
Ichiye  [16]  « We  shall  discuss  his  work  later  in  the  report. 

It  has  been  generally  agreed  upon  by  oceanographers 
that  the  type  of  phenomena  we  wish  to  consider  can  be  adequately 
described  by  the  dynamics  of  the  problem  alone,  the  temperature 
effects  being  included  by  way  of  an  assumed  semi-empirical  den- 
sity distribution.  At  the  Woods  Hole  Oceanographic  Institute, 
experiments  with  a model  parabolic  ocean  basin  verify  the  above 
conjecture.  Hence,  in  the  subsequent  analysis,  we  shall  neglect 
direct  temperature  dependency  in  the  treatment  of  the  problem 
and  shall  Include  only  the  effects  of  wind  and  gravitation. 

A large  part  of  our  report  is  concerned  with  the  formu- 
lation of  the  problem  and  the  assumptions  made  to  reduce  the 
general  problem  to  one  which  can  be  attacked  mathematically.  In 
the  past  a discussion  of  such  assumptions  has  often  been  vague. 
It  was  felt  therefore  that  an  explicit  and  detailed  analysis  of 
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the  simplifications  involved  in  the  formulation  of  the  problem 
might  be  welcomed  by  workers  in  this  field  and  that  it  might 
help  to  clear  up  any  existing  misconceptions  concerning  the 
validity  of  some  of  the  assumptions. 

2.  Discussion  of  Results.  At  this  point  we  shall  discuss, 
without  resorting  to  mathematical  detail,  the  basic  assumptions, 
the  results,  and  the  conclusions  of  the  present  investigation. 

In  this  manner  we  hope  to  convey  a more  integrated  picture  of 
the  physics  of  the  problem. 

Mathematically,  the  motion  which  we  want  to  study  can 
be  defined  by  the  Navier-Stokes  equations  of  motion  with  the 
viscous  terms  replaced  by  terms  arising  from  a macroscopic  vis- 
cosity, viz, , an  eddy  viscosity.  The  complete  non-linear  equa- 
tions are  too  difficult  to  solve,  however,  so  that  we  are  forced 
to  make  a number  of  simplifying  assumptions  which  we  shall  list 
below. 

1.  The  fluid  is  assumed  to  be  incompressible,  but  it  may 
be  Inhomogeneous. 

2.  The  equations  on  a rotating  sphere  are  approximated  by 
equations  in  a rectangular  Cartesian  system.  The  effect  of  the 
sphericity  of  the  earth  is  retained  by  allowing  the  Coriolis 
parameter  to  depend  on  the  latitude.  Since  ye  shall  consider  a 
rectangular  ocean  in  the  Cartesian  system,  a few  remarks  must 
be  made  concerning  the  region  of  the  sphere  onto  which  the  rec- 
tangle is  mapped.  The  constant  east-v;est  distance  of  the  rec- 
tangle is  preserved  in  the  mapping  of  the  rectangle  onto  the 
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spheret  Such  a mapping  is  not  conformal  since  angles  between 
lines  are  not  preserved*  The  region  under  consideration  must 
be  well  removed  from  the  north  pole* 

3*  The  vertical  acceleration  terms  and  the  viscous  terms 
are  neglected  in  the  equation  of  vertical  motion  so  that,  in 
effect,  hydrostatic  pressure  is  assumed,  i.e.,  p = gj^pdz,  where 
T]  is  the  free  surface  height  and  p = 0 at  z = t].  The  density 
p may,  of  course,  be  a function  of  the  space  coordinates*  In 
Appendix  3 it  is  shown  that  for  the  problem  which  is  independ- 
ent of  time,  the  hydrostatic  pressure  assumption  is  necessary 
only  in  the  depths  where  there  is  no  motion  if  one  desires  a 
solution  for  the  components  of  the  mass  transport  only,  If  it 
is  necessary  to  find  the  shape  of  the  free  surface,  however,  or 
if  the  non-steady  problem  is  considered,  this  assumption  or  some 
analogous  one  must  be  made* 

4*  As  stated  in  the  introduction,  the  thermodynamic  effect 
are  accounted  for  only  empirically  by  stipulating  a density  dis- 
tribution* We  assume  p = p[z  - T(x,y,t)]  where  the  function  p 
of  the  variable  (z  - T)  can  be  prescribed  to  fit  observational 
data*  This  functional  form  for  p makes  the  curves  of  constant 
density  parallel* 

5*  The  equations  of  motion  are  integrated  over  the  verti- 
cal coordinate,  z* 

In  order  to  perform  this  integration  it  is  necessary 
that  we  specify  the  density  distribution  since  p appears  in 
some  of  the  integrands*  V/e  consider  two  cases, 

(1)  The  surface  z = T separates  two  layers  of  constant 
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density.  The  equations  of  motion  in  each  layer  are  then  inte- 
grated over  the  depths  of  the  respective  layers  and  the  non- 
linear terms  are  neglected.  We  also  neglect  shear  forces  at 
the  bottom  of  the  lower  layer  and  at  the  interface.  No  assump- 
tion  is  made  concerning  the  vertical  distribution  of  velocity  , 
but  instead,  we  hope  to  solve  for  the  integrated  velocities 
(i,e. , the  transports)  in  each  layer.  This  case  is  referred  to 
as  the  two-layer  problem.  Unfortunate ly^  it  is  much  too  diffi- 
cult to  handle  analytically,  and  consequently  we  must  consider 
a second  problem. 

(ii)  The  manner  of  performing  the  integration  in  this  case 
will  lead  to  a considerably  simplified  problem  which  allows  us 
to  stipulate  a more  general  density  distribution  than  that  in 
(i).  The  density  is  specified  as  a continuous  function  of  depth 
and  the  ocean  is  divided  into  three  layers,  A layer  of  constant 
density,  pQ,  lies  above  the  surface  z = T(x,y,t).  From  z = T 
down  to  z = T - d (d  is  constant)  the  density  increases  linearly 
with  depth  from  Pq  to  the  value  p Below  z = T - d,  the  den- 
sity has  the  constant  value,  p 

We  assume  that  there  is  a depth  z = - h(x,y,t)  belov; 
which  the  velocities  may  be  considered  negligible  (in  some 
suitably  defined  sense).  The  pressure  gradients  will  then  also 
be  negligible  below  z = - h.  As  a consequence  of  this  assumption 
and  the  previous  assumption  of  hydrostatic  pressure,  a relation- 
ship exists  between  the  surface  z = T and  the  free  surface 


* Compare  this  xvith  case  (iiTl. 
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z = T]  , Viz.  , T = - t)  - C (where  Ap  = P **  P o T = -C 

when  T]  = 0).  Thus,  If  the  velqclties  are  negligible  In  the 
depths  of  the  ocean,  the  thermocline  must  respond  immediately 
to  a change  in  the  shape  of  the  free  surface  in  order  to  main- 
tain negligible  pressure  gradients  at  these  depths. 

The  three  assumptions,  (a)  hydrostatic  pressure, 

(b)  negligible  velocities  in  the  ocean  depths,  and  (c)  con- 
stant density  below  the  thermocline,  are  crucial  for  the  present 
case.  It  is,  of  course,  possible  that  any  one  or  a combination 
of  these  three  assumptions  may  be  incorrect.  If  this  be  the 
case,  then  the  thermocline  need  not  respond  to  the  free  surface 
immediately.  The  frequency  of  the  wind  variation  which  we  shall 
consider  later  in  our  development  will  be  small  so  that  assump- 
tions (a)  and  (b)  seem  plausible.  . Thus  the  only  motion  exist- 
ing below  the  thermocline  is  caused  by  vertical  shear  and  this 
motion  decays  exponentially  with  increasing  depth  according  to 
Ekman  [ l]  • 

The  equations  of  motion  are  then  integrated  from  the 
depth  z = - h to  the  free  siirface  z = q • This  problem  will  be 
called  the  one-layer  problem  because  of  the  single  integration. 
The  depth,  z = - h,  does  not  appear  explicitly  in  the  Integrated 
equations. 

In  both  cases,  the  effect  of  the  wind  is  represented 
by  the  shear  stress  at  the  ocean  surface  and  appears  in  the 
evaluation  of  the  vertical  viscous  terms  at  the  upper  limit  of 
integration  (free  surface). 

An  additional  difference  between  the  two  problems  is 
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that  the  two-layer  problem  specifically  restricts  the  fluid  of 
the  top  layer  to  remain  in  the  top  layer  and  the  fluid  in  the 
lower  layer  to  remain  In  the  lower  layer#  The  one-layer  prob- 
lem has  no  such  restriction  and  an  interchange  of  fluid  may 
result.  However,  because  of  the  integration  we  have  no  inform- 
ation concerning  this  vertical  motion. 

6.  The  non-linear  terms  in  the  equations  of  horizontal 
motion  are  neglected,  A plausibility  argument  for  this  assump- 
tion, based  on  the  results  of  [7]>  is  presented  in  Appendix  2. 
However,  our  results  must  now  be  considered  tentative,  since 
the  case  presented  in  the  appendix  for  the  neglect  of  the  non- 
linear terms  is  a plausibility  argument  and  not  a justification# 
The  primary  motive  for  neglecting  the  non-linear  terms  is  our 
inability  to  cope  with  them  analytically# 

7.  The  Coriolis  parameter  is  linearized.  In  effect, 
this  is  comparable  to  linearizing  the  sine  of  an  angle  when  the 
angle  varies  between  l5°  and  60°, 

With  the  above  assumptions  and  simplifications  we  are 
in  a position  to  attempt  a solution  of  the  non- steady  problem. 
The  ocean  is  chosen  to  be  rectangular  with  vertical  walls  as 
boundaries  on  the  east  and  west.  Because  of  the  presence  of 
viscosity,  the  boundary  conditions  on  these  walls  are  that  the 
velocities  vanish#  The  boundaries  on  the  north  and  south  are 
water  boundaries# 

The  wind-stress  is  v/ritten  as 

= - (VJ*  +r’sinwt)cos  ny 
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where  W*,P  *,  w , and  n are  constants  and  (Fig.  1)  is  the 
east-west  component  of  the  stress.  The  above  form  for  the  wind- 
stress  may  be  considered  as  the  general  term  of  a Fourier  series 
expansion  so  that  the  wind-stress  may  be  generalized  for  the 
linear  problem.  However,  for  our  numerical  example,  we  have 
chosen  to  to  give  a period  of  one  year  and  n as  2u/s  where  s is 
the  north-south  length  of  the  ocean  (0  ^ y s).  The  wind- 
stress  component  Ty  is  assumed  Identically  zero.  Since  the  wlnd- 
stress  is  prescribed  in  such  a manner  that  its  y derivative 
vanishes  at  y = 0,s,  it  appears  reasonable  to  demand  that  these 
boundaries  be  streamlines  and  that  the  normal  derivatives  of 
the  velocities  vanish  there. 

The  one-layer  problem  is  solved  by  the  following  proce- 
dure. The  equations  are  non-dimensionalized.  The  non-dimen- 
sional velocities  and  free  surface  height  are  expanded  in  per- 
turbation series  with  the  non-dimensional  time  parameter  as  the 
perturbation  parameter.  Each  resulting  set  of  equations  is 
then  solved  by  application  o'f  the  boundary  layer  technique* 

The  conditions  for  the  validity  of  the  expansion  restrict 
the  time  variation  to  a maximum  frequency  of  seasonal  oscilla- 
tion, In  the  numerical  example,  yearly  frequency  is  assumed 
and  the  perturbation  terms  of  second-order  and  higher  are 
neglected.  The  error  Involved  in  neglecting  the  second-order 
term  as  compared  to  the  zero-order  term  is  about  and  it  is 
about  20%  as  compared  to  the  first-order  term.  The  remaining 
physical  parameters  are  given  values  which  correspond  roughly  to 
those  of  the  North  Atlantic  Ocean, 
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The  following  discussion  will  he  based  on  the  non- 
dimensional  quantities  defined  in  the  body  of  the  report.  When- 
ever dimensional  quantities  are  mentioned,  we  shall  include  the 
dimensions. 

The  graph  of  the  north-south  component,  V,  of  the  mass 
transport  vs.  the  east-west  coordinate  x’  near  x’  = 0,  the 
western  shore,  is  shown  in  Fig.  2 for  the  value  y*  = 0.25,  i.  e,  , 
where  the  Gulf  Stream  is  most  pronounced.  The  Gulf  Stream  re- 
gion is  the  region  of  large  positive  V.  The  region  of  negative 
V adjacent  to  the  Gulf  Stream  corresponds  to  the  offshore 
counter-current. 

The  Gulf  Stream  responds  to  the  wind  in  such  a manner 
that  the  mass  transport  and  the  wind  are  in  phase  whenever  the 
latter  takes  on  its  maximum  or  minimum  value.  At  all  other 
times  the  mass  tramsport  lags  behind  the  wind  with  the  greatest 
lag  occurring  when  the  wind  reaches  its  steady  position  . At 
this  time  the  mass  transport  is  about  9 days  away  from  its 
steady  value  . The  length  of  this  InVeryal,  i.e»  , nine  day:s, 
is  independent  of  the  frequency  for  slowly  varying  winds. 

The  wind  (see  Fig.  1)  and  the  mass  transport  attain 
their  maximum  values  at  t = n/2.  The  mass  transport  now  has  a 
magnitude  of  (1  + r»/W’)  times  its  steady  value.  Thus,  within 
the  accuracy  of  the  present  method  of  solution,  the  time  at 
which  maximum  transport  occurs  and  the  magnitude  of  the  maximum 

* We  shall  refer  to  the  "steady  position"  whenever  the  time- 
dependent  contribution  of  the  wind  is  zero. 

**  l.e.  , the  value  due  to  its  response  to  a steady  wind 
= - W ‘ cos  ny. 
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transport  are  independent  of  the  frequency.  The  magnitude  of 
the  out-of-phase  effect  (the  second  term  in  the  perturbation 
series)  which  is  largest  when  the  wind  has  its  steady  value, 
is  proportional  to  the  frequency. 

The  time  variation  of  the  wind  affects  the  Gulf  Stream 
only  by  changing  the  mass  transport  through  the  Stream.  It 
does  not  change  the  Stream* s position. 

As  can  be  seen  from  Fig,  2,  the  relative  importance  of 
the  out-of-phase  effect  is  greatest  in  the  counter-current. 

Figure  3 is  a graph  of  the  north-south  mass  transport 
component  near  the  eastern  boundary  of  the  rectangular  ocean 
at  the  latitude  y*  = 0,25,  The  accompanying  out-of-phase  effect 
is  shown  at  its  maximum  in  the  figure,  V is  negative  on  the 
eastern  coast,  i,e.,  the  mass  transport  is  toward  the  south. 

Figures  4,  5,  and  6 show  the  contour  lines  of  the  free 
surface  in  the  southern  half  of  the  ocean  for  various  times* 

With  the  values  of  the  contour  lines  multiplied  by  -200  the 
three  figures  represent  the  contour  lines  of  the  thermocline. 
Qualitatively,  the  results  agree  fairly  v/ell  with  observation 
though  some  of  the  natural  features  are  missing.  It  seems 
likely,  however,  that  most  missing  features  result  from  local 
effects  v/hich  we  have  not  taken  into  account. 

Because  of  the  lengthy  computations  Involved,  we  have 
calculated  numerical  results  for  only  one  set  of  values  of  the 
parameters.  It  can  be  seen  from  the  analytical  results  that  if 
the  average  depth  of  the  top  layer  be  changed,  the  values  for 
the  deflection  of  the  free  surface  and  the  out-of-phase 
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velocities  will  change,  Specifically,  if  the  depth  is  decreased, 
the  free  surface  deflection  is  increased  and  all  out-of-phase 
quantities  are  increased. 

The  above  results  appear  to  invalidate  the  solution  of 
the  problem  as  obtained  by  Ichiye  [l6]  • Ichiye  neglected  the 
contribution  of  the  non-steady  term  in  the  Integrated  continuity 
equation.  However,  with  the  values  of  the  parameters  used  in 
Section  the  magnitude  of  this  term  in  the  Interior  of  the 
ocean  is  as  much  as  ten  times  that  of  the  remaining  non-steady 
terms  which  were  retained  in  Ichiye 's  analysis* 

Vie  have  computed  the  mass  transport  through  the  Gulf 
Stream  for  the  one-layer  steady  problem.  With  the  given  wind 
distribution  our  result  is  26,6  x 10°  metric  tons  per  second. 

This  value  is  about  three-fourths  of  Munk’s  value  [^]  and  about 
one-third  of  the  observed  value,  Munk  used  an  empirical  east- 
west  wind  distribution. 

The  two-layer  steady  problem  is  solved  in  Section  ? 
where  it  is  shown  that  the  mass  transport  streamline  pattern  is 
the  same  as  in  the  one-layer  problem.  This  is  to  be  expected 
since,  for  the  steady  case,  the  same  assumptions  are  made  regard- 
ing negligible  velocities  below  the  thermocllne.  Thus,  the 
height  of  the  thermocllne  is  shown  to  be  proportional**  to  the 
free  surface  deflection.  Since  the  free  surface  height  is  deter- 
mined largely  by  the  thickness  of  the  top  layer,  the  thermocllne 

* In  [l6J  the  term  corresponding  to  W»  in  the  present  paper 
was  assumed  to  be  Identically  zero,  1, e, , the  wind  had  a 
zero  mean  value, 

**  The  factor  of  proportionality  is  the  reciprocal  of  the 
density  difference. 
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variation  depends  on  the  choice  of  the  two  parameters,  density 
difference  and  thickness  of  top  layer* 

By  varying  the  two  parameters  we  can  get  good  qualita- 
tive agreement  with  observations  of  the  shape  of  the  thermocline. 
In  Fig,  9,  a cross-section  of  the  computed  thermocline  is  shown 
for  four  pairs  of  values  of  the  parameters.  Because  of  the 
rather  vague  definition  of  the  actual  thermocline,  we  cannot 
state  specifically  the  extent  of  quantitative  agreement  between 
our  computed  results  and  the  observed  values.  Consider, 
however,  the  curve  in  Fig,  9 with  a depth  of  the  top  layer  of 
200  meters  and  a density  difference  of  0,0025#  For  that  curve 
the  results  disagree  by  a factor  of  three  when  compared  to  some 
of  the  measurements  of  the  thermocline  off  Chesapeake  Bay  [lO], 
The  two-layer  non-steady  problem  constitutes  an  attempt 
to  drop  the  assumption  made  in  the  one-layer  problem  that  the 
velocities  vanish  at  some  great  depth.  As  a consequence  the 
problem  becomes  much  more  complicated  and  it  is  necessary  to 
introduce  some  other  simplifying  assumptions,  viz. , to  neglect 
the  shear  forces  at  the  bottom  and  at  the  thermocline.  This 
may  have  far-reaching  effects.  These  simplifications  notwith- 
standing, we  were  unable  to  obtain  a solution,  A brief  descrip- 
tion of  our  attempts  at  such  a solution  follows. 

First,  the  equations  are  non-dimensionalized  as  in  the 
one-layer  case.  The  integrated  continuity  equation  for  the  top 
layer  now  contains  the  time  derivative  of  the  magnitude  of  the 
deviation  of  the  thermocline  from  its  equilibrium  position. 

Since  this  term  is  very  large,  the  perturbation  method  used  in 
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the  one-layer  problem  Is  restricted  to  a range  of  frequency 
values  corresponding  to  less  than  one  oscillation  every  hundred 
years*  Since  these  results  are  not  physically  interesting  no 
numerical  results  were  computed* 

A second  method  of  attack  is  then  attempted#  The  wind- 
stress  term  is  first  divided  into  its  steady  and  non-steady  parts 
and  the  two  problems  are  treated  separately  without  resorting 
to  a perturbation  in  the  time  parameter*  This  method  had  been 
attempted  for  the  one-layer  problem  with  no  success*  In  the 
present  case,  however,  it  was  hoped  that  the  new  parameter  in- 
volving the  density  difference  could  be  used  to  advantage.  Un- 
fortunately, an  analytic  solution  still  appears  to  be  quite 
hopeless* 

The  one  interesting  fact  which  seems  to  emerge  from  the 
attempts  at  the  solution  of  our  idealized,  two-layer,  non-steady 
problem  concerns  the  magnitude  of  the  lower  layer  transport*  We 
must  recall  that,  in  the  case  treated,  the  solution  is  restricted 
to  the  frequency  range  for  which  the  thermocllne  responds  to  the 
variation  of  the  top  surface  in  a quasi-steady  manner;  i. e. , as 
a result  of  any  change  in  the  free  surface,  the  thermocllne 
assumes  the  same  shape  as  it  would  for  a steady  problem  with  the 
given  free  surface,  except  for  a small  out-of -phase  correction# 

In  this  case,  the  mass  transport  in  the  lower  layer,  excluding 
whatever  transport  may  be  caused  by  shear  at  the  Interface,  is 
of  the  same  order  of  magnitude  as  that  portion  of  the  transport 
in  the  upper  layer  which  is  out  of  -phase  with  the  wind#  For  a 
higher  frequency  this  result  does  not  necessarily  hold  true* 
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A final  v/ord  should  be  said  about  the  lack  of  quantita- 
tive agreement  between  our  computed  results  and  observation. 

The  factor  of  three  is  not  surprising  when  one  considers  the 
very  idealized  model  which  we  have  assuiiied,  A number  of  more 
realistic  assumptions  may  certainly  affect  our  quantitative 
results  by  such  a factor.  The  inclusion  of  the  non-linear  terms, 
a better  representation  of  the  wind  effects  on  the  water,  a 
more  natural  topography,  and  a non-constant  eddy  viscosity  may 
well  alter  the  quantitative  results  and  bring  them  into  closer 
agreement  with  reality. 

3*  Formulation  of  the  Problem.  It  is  our  aim  to  derive 
expressions  for  the  velocity  and  the  pressure  satisfying  the 
three  equations  of  motion  on  a rotating  sphere 

9q  T T 

■^  + q • + 2S2xq  + Qx(Qxr)  = - -*  Vp  + F + -(V  * A,-V)q 

at  — — — p — p — 

the  continuity  equation 

V • q = 0 

and  the  boundary  condition  that  q = 0 on  a land-water  boundary* 
Here,  ^ 

q = (u,v,v/)  denotes  the  velocity  vector  relative  to  a 
“ coordinate  system  rotating  with  the  sphere, 

Q,  denotes  the  angular  velocity  vector  representing  the 
’’  earth's  rotation, 

p denotes  the  pressure, 

p denotes  the  density, 

F denotes  the  external  forces  per  unit  mass  (in  our  case, 
” gravitation), 

* u,v,w  are  spherical  components  of  velocity  along  the  direc- 
tions of  the  radius,  the  meridians,  and  the  parallels  of 
latitude  respectively* 
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(7  • A.V)q  represents  the  eddy  viscosity  term  (discussed 
below). 

Let  us  consider  the  expression  for  the  eddy  viscosity 
term  in  a rectangular  coordinate  system,  this  being  the  system 
in  which  we  shall  later  write  our  equations. 

We  define  the  operator  (7  • A^7)  as  follows; 


where  A^^,  A2,  A^,  may  depend  on  the  space  coordinates.  These 
three  quantities  (the  coefficients  of  the  lateral  and  vertical 
eddy  viscosity)  have  been  measured  and  are  knovm  to  vary  through- 
out the  ocean..  The  definition  of  the  viscous  coefficients  and 
our  knowledge  of  their  magnitudes,  however,  are  rather  vague* 

In  view  of  this,  and  because  of  subsequent  analytical  simplifi- 
cations, we  assume  that  the  lateral  kinematic  eddy  viscosity 
coefficients  are  constant  and  equal,  so  that 


Where  A is  now  a kinematic  eddy  viscosity  and  is  constant.  No 
simplification  will  be  made  concerning  A^* 

Our  continuity  equation  is  valid  for  an  Incompressible 
fluid.  In  the  steady  problem  the  density  may  be  more  general 
and  we  have  simply  7 • (pq)  = 0,-  In  the  non-steady  problem,  the 
assumption  of  incompressibility  is  imposed  but  the  fluid  may  be 
homogeneous. 

We  shall  want  to  make  use  of  [7]  regarding  the  effect 


of  the  non-linear  terms.  Because  the  results  in  [7]  are  discussed 
in  terms  of  rectangular  coordinates  and  because  the  use  of 
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rectangular  coordinates  considerably  simplifies  the  analysis, 


we  shall  first  transform  the  equations  of  motion  from  spherical 
to  rectangular  coordinates  in  such  a manner  that  the  equilibrium 
free  surface  which  establishes  itself  in  the  spherical  system 
as  a result  of  gravity  and  centripetal  acceleration  corresponds 
to  the  x-y  plane  of  the  rectangular  system.  The  apparent  gravi- 
tational force,  i.e. , the  force  which  is  the  resultant  of  true 
gravity  and  centripetal  acceleration,  acts  in  a direction  normal 
to  this  equilibrium  surface* 

In  Appendix  1,  it  is  shown  that  our  original  equations 
reduce  to 


+ u»  + v«  In!  - 2QV*  sin(I)  s - 1 M+4(V»A,V)u‘  (1) 

8t  ox  3y  R |3  ax  J 1 


+ u‘  + v»  ^ + 2Qu«  sin(|)  = - i iPU(VA.V)v»  (2) 

at  3x  ay  R p ay  > i 


- 1 M 


p 9z 


g 


(3) 


where 


ml  + = 0 

ax  ay  az 


(4) 


x, u’  denote  the  east-west  coordinate  and  velocity 

respectively  (x  is  positive  eastward), 

y, v'  denote  the  north-south  coordinate  and  velocity 

respectively  (y  is  positive  northward), 

z, w'  denote  the  vertical  coordinate  and  velocity 

respectively  (z  is  positive  upward), 


R is  the  mean  radius  of  the  earth. 


g is  the  apparent  gravitational  acceleration  on 
the  earth’s  surface, 

2Qsln(^)  is  the  radial  com.ponent  of  the  angular  velocity 
^ vector  of  the  earth. 
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The  rectangular  coordinate  system  is  oriented  with  the 
origin  in  the  southwest  corner  of  the  ocean  and  with  the  equili- 
brltani  surface  in  the  x-y  plane. 

A nvunber  of  assumptions  were  made  in  the  reduction  of 
the  four  equations  valid  on  a spherical  earth  to  the  four  equa- 
tions given  above^  These  assumptions  are  listed  here  for  the 
convenience  of  the  reader  who  does  not  wish  to  go  through  the 
detail  in  Appendix  I4 

(1)  In  the  radial  component  of  the  equations  of  motion^ 
the  acceleration  terms  and  the  viscous  terms  are 
neglected  in  comparison  to  g,  the  gravitational  accel- 
oration.  In  essence,  we  assume  hydrostatic  pressure  • 

(2)  All  terms  involving  radial  velocity  are  neglected  in 
the  remaining  two  equations  of  motion  on  the  supposi- 
tion that  the  radial  velocity  is  very  small  compared 
to  the  lateral  velocities. 

(3)  The  variation  of  the  radial  distance,  r,  over  the 
depth  of  the  ocean  is  neglected  and  we  write  r ^ R, 
the  mean  radius  of  the  earth* 

(Actually^  tho  radial  distance  varios  by  about  1/1000 
of  its  total  length. ) 

(*f)  Terms  which  are  divided  by  R are  neglected  in  compar- 
ison with  all  other  terms. 

(5)  The  region  considered  must  not  lie  close  to  the  north 
pole  since  some  terms  which  have  been  neglected 


♦ In  Appendix  3j  this  assumption  is  discussed  in  more  detail# 
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previously  become  infinite  at  the  pole.  In  our  prob- 
lem the  ocean  is  confined  to  a region  lying  south  of 
latitude  70°. 

(6)  An  appropriate  interpretation  of  the  results  as  applied 
to  the  spherical  earth  must  be  made,  keeping  in  mind 
that  the  boundaries  have  been  distorted.  If  we  con- 
sider a rectangular  ocean  in  the  plane,  the  appropri- 
ate mapping  onto  the  sphere  v/ould  preserve  the  con- 
stant east-west  length.  Such  a mapping  is  not  conform- 
al siiice  angles  are  not  preserved,  (In  the  case  .of  a 
Mercator  projection,  on  the  other  hand,  angles  are 
preserved,  bub  the  east-west  distance  is  distorted.) 

Let  us  consider  the  simplified  equation  of  vertical 
motion  (3)*  In  integrated  form,  this  equation  is 

p = g pdz  (3»a) 

Jz 

where  q measures  the  deflection  of  the  free  surface  from  its 
equilibrium  position  and  the  scale  of  p is  chosen  in  such  a 
manner  that  p = 0 on  z = q.  Wov/,  the  density  is  a function  of 
temperature  and  salinity.  In  our  treatment  of  the  problem,  how- 
ever, v/e  wish  to  avoid  the  analytical  difficulties  introduced 
by  Including,  explicitly,  the  energy  eqTiation  and  an  equation  of 
state,  \/e  propose  irste.ad  to  account  for  the  thermodynamics  of 
the  problem  empirically  by  prescribing  a density  distribution 
which  roughly  conforms  to  observation’^.  In  particular,  we 

* In  Appendix  3 it  is  st ovm 'that  a specification  of  the  density 
distribution  and  the  tssuraption  of  hydrostatic  pressure  are 
not  necessary  for  the  steady  problem# 
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choose  p = p[z  - T(x,y,,t)],  where  the  function  p of  the  variable 
(z  - T)  can  be  prescribed  to  fit  observational  data.  We  observe 
that  this  functional  form  for  p makes  the  curves  of  constant 
density  parallel  to  each  other. 

A complete  analysis  for  the  unknown  quantities  as  func- 
tions of  the  four  independent  variables  x,y,z,t  is  exceedingly 
difficult  and  we  are  forced  to  eliminate  one  variable  by  inte- 
grating our  equations  over  the  vertical  coordinate,  z,  and  then 
solving  for  suitably  defined  integrated  quantities#  In  so  doing, 
we  lose  information  concerning  the  dependence  of  the  unknowns 
on  z.  Since  we  are  primarily  concerned  with  general  oceanic 
circulation  and  mass  transport,  however,  and  since  the  integra- 
tion leads  to  a considerable  reduction  in  dif faulty,  the  advan- 
tages gained  more  than  balance  the  loss  of  information  Involved. 

Actually  we  cannot  afford  a complete  loss  of  information 
concerning  the  vertical  dependence  of  velocity.  This  will  become 
apparent  shortly. 

The  general  density  distribution  must  be  specialized  in 
order  to  permit  Integration  of  the  equations  over  the  vertical 
coordinate.  Two  cases  v/111  be  considered. 

First 5 let  T be  a surface  which  separates  two  layers  of 
constant  density  so  that 


p[z  - T(x,y,t)]  = p2  for  z > T(x,y,t) 

and 

p[z  - T(x,y,t)]  = P2  for  z < T(x,y,t). 

For  this  problem  it  is  convenient  to  choose  the  coordi- 
nate system  with  the  xy-planes  parallel  to  the  undisturbed 
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equilibrium  surface  and  \]±th  the  plane  z = 0 at  the  bottom  of 
the  ocean,  the  bottom  being  assumed  plane  in  this  problem,  A 
layer  of  constant  density  P2  extends  from  the  bottom  of  the 
ocean  to  the  height  z = D2  + ^2  where  the  constant  D2  is  the 
average  height  of  the  lower  layer  and  r\2  is  the  height  of  the 
disturbed  surface  of  this  layer  measured  from  the  plane  z = D2, 
A layer  of  constant  density  pj^  extends  from  the  height 
z = D2  + TI2  surface  z = where  is  the 

distance  from  z = 0 of  the  undisturbed  equilibrium  surface  of 
the  upper  layer  and  r)]_  is  the  height  of  the  disturbed  free  sur- 
face of  the  upper  layer  measured  from  z = Dp 
Then  equation  (3«a)  becomes 


Pt  = gPi(  1i  + - z)  for  the  upper  layen 

^ (3.b) 

P2  = 6Pi(  ^ + ^1  " ^2  " ^2^  SP2(^2  + ^2  “ lower  layes 

(3.0) 

If  we  denote  all  quantities  in  the  upper  and  lower  layers 
by  subscripts  1 and  2,  respectively,  the  equations  (1),  (2)  and 
(If),  with  expressions  (3«b)  and  (3.0)  substituted  for  the  pres- 
sure in  the  upper  and  lower  layers,  respectively,  become 


•5^  * “1  Ir  ^ ^ ® ^ ^ 


8v 

at 
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(9) 

(10) 


where  a = b =(P2  “ P^/P2  =^p/p2* 

The  problem  defined  by  equations  (5)  - (10)  with  appro- 
priate boundary  conditions  is  quite  general  in  that  no  assumption 
has  been  made  concerning  the  vertical  distribution  of  velocity. 

As  we  shall  see  later,  when  the  equations  are  integrated  over  z 
and  linearized,  the  simplified  problem  is  still  too  difficult 
to  solve.  For  this  reason  we  formulate  a second  problem  v/hich 
allows  a more  general  density  distribution  but  which  is  more 
restricted  in  other  respects* 

In  this  problem  we  retain,  for  the  time  being, the  gen- 
eral form  p = p [z  - T(x,y,t)].  Then  the  pressure  terms  in 

die 

equations  (1)  and  (2)  are 


I 


I ^ S 

p ax  p 


p ay  p 


(11.  a) 


(ll.b) 


* For  the  present  problem  the  plane  z = 0 lies  on  the  undisturb- 
ed equilibrium  free  surface. 
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where  Pq  = p[fi  - T(x,y,t)  ],  the  density  at  the  free  surface. 

If  these  terms  be  substituted  into  (1)  and  (2),  we  have 


^ + u»  + v'  - 2Qv‘  sin  (?)  = 
ot  0x  ay  R 


u 


£ dz  - ^ p^  +i(V  • A.V)u‘  (12) 

X p 6x  0 f ^ 


9p 

8x 


^ + u»  ^ + V'  ^ + 2Tm»  sin  (|)  = 
at  6x  ay  R 


g 


3p 


p ^z  9y 


dz  - i £n  p^  +jp(  V • A.V)v».  (13) 

p ay  o f 1 


As  stated  previously,  the  problem  will  be  simplified  by 
Integrating  the  equations  over  the  vertical  coordinate,  z. 

Let  us  first  consider  the  problem  defined  by  the  equa- 
tions (k) , (12),  (13).  We  assume  that  there  is  a depth 
z = - h(x,y,t)  below  which  the  velocities  may  be  considered 
negligible  (in  some  suitably  defined  sense),  and  we  integrate 
from  z = - h up  to  the  free  surface.  The  depth  z = - h(x,y,t) 
may,  of  course,  vary  from  point  to  point  in  the  ocean.  Since 
the  velocities  are  negligibly  small  below  z = - h,  the  horizontal 
pressure  gradients  must  also  be  negligibly  small  and  we  may 
therefore  write 


1 ^ 

p 3x 


= 0. 


z=-h 


I ^ 
p ay 


= 0, 

z=-h 


(14) 


V/e  must  now  specialize  the  general  form  of  the  density 
distribution  because  an  integration  involving  p will  actually 


* This  assumption  is  the  fundamental  difference  between  the 
two  problems  considered. 
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have  to  be  carried  out. 

Define  p = p[z  - T(x,y,t)]  in  such  a way  that 

p = Pq,  a constant  for  > 2 > ^ 
p = [1  + c(T  - z)]pq  for  T 2 z 2 constants) 

P = p_j^  = (1  + cd)pQ  for  T - d > z . 


r 

r (15) 


With  this  definition  the  density  is  a continuous  function  of 
depth  and  the  ocean  is  divided  into  three  distinct  layers,  A 
layer  of  constant  density,  p^,  lies  above  a region  in  which  the 
density  increases  linearly  with  depth  from  p^  to  the  value  p_j^. 
Finally,  at  the  bottom,  there  is  a layer  of  constant  density, 
p_^.  This  prescribed  distribution  agrees  well  with  the  observed 
density  distribution. 

If  p,  as  given  by  (1^),  be  substituted  into  equation 
(14),  we  find  that* 

^ = _£olH  1T  = -£.o£II  (16) 

3x  ” Ap  9x’  3y  Ap  3y  '■ 

where  Ap  = P_h  ” Po* 

If  we  integrate  equations  (I6),  v/e  obtain 

T = - n - c (17) 

Ap 

where  z = - C is  the  constant  depth  of  T when  q = 0,  Physically, 
z = - C is  an  average  depth  of  the  top  layer  or  the  depth  of  T 
when  the  ocean  surface  is  undisturbed  (i.e, , in  the  absence  of 
winds).  These  tv/o  quantities  are,  of  course,  identical. 


* The  algebraic  manipulation  is  given  in  Appendix  4(a), 
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Let  us  next  integrate  equations  (12)  and  (13)  from 
h to  z = T]t  The  pressure  terms  become 


(18,  a) 


1 |P  dz  = - gD 
P 9y 


Sn  . g „ 2a 

9y  ^ Ap  ^ 9y 


where  D = C + d/2,  and  the  complete  equations  are 


(I8,b) 


9U  +-T  9ul  (i2  + p v«  MI  dz  - 2T^V  sin  (|) 

at  ^ ax  ^ dy  R 

J-h  i^-h 


= - gD  Me.  - g ^ n Me.  +AAU  + (A 
ax  Ap  ^ ax 


3 9z 


(19) 


- - r ^ rn 

av+p  u»MIdz+p  v»  MI  dz  + 2S2U  sin  (|) 

^ 9x  9y  R 

J -h  ^-h 

= . gD  . gl=h  |3£  + + (A,  |^)j  ’’ 

9y  Ap  ‘ ay  3 az 

where  _ f M — f 

U = pu’dz,  V = ^v'dz, 

J -h  -h 

^ is  a constant, average  density, 

and  \(x,y,z,t)  j = \(x,y,q,t)  - \(x,y,-h,t). 

The  non-linear  terms,  u'(x,y,  T],t)9q/at,  etc.,  from 


(20) 


* See 


ihe  details. 


Since  the  viscous  terms  are,  in  any  case,  only  approximations 
to  the  actual  shear  stresses,  we  have  made  the  further  approx- 
imation 

f ^ (At  .2iil)dz  -if  9 / A au *1  — 1a  au*  I ^ 

J -h  p 82  3 az  J5-  J .h  ^^*3  arldz  r 3 5^  |_^  • 
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the  interchange  of  integrals  and  derivatives  of  the  velocity- 
terms  have  been  neglected,  ¥e  have  defined  U and  V as  mass 
transport  components  rather  than  as  volume  transport  components 
(by  simply  including  an  average  density  in  the  definition)  be- 
cause we  want  to  compare  some  of  our  quantitative  results  with 
observations  and  with  the  results  of  Munk^  both  of  which  are 
given  in  terms  of  mass  transport. 

The  terras  liil  1 ^ and  Ao  9v  * I ^ must  give  the  wind- 

stress  terms  since  they  represent  the  shear  stress  evaluated  at 
the  upper  limits  (the  shear  stress  terms  at  z = - h are  negli- 
gible since  -h  was  chosen  as  the  depth  v/here  the  motion  becomes 
negligible).  Thus 

A-;  -231-  I = = X component  of  wind  stress 

b 8z 

Ao  I ^ = y component  of  wind  stress# 

■j  oz  y 

In  the  equation  of  continuity  we  shall  want  to  make  use 
of  the  kinematic  free  surface  condition  [9] 


^ [z  - n(x,y,t)]  =0  at  z = q • 

When  expanded,  this  equation  reads 

= ^ + u»p  ^ + v»|  ^ 4^ 

'at  'ax  ' ^ 

where  v7 ' P etc,  denotes  the  value  of  w *(x,y,z,t)  at  z = q , 
Integration  of  the  continuity  equation  (4)  yields 


pv'pAq  + pw’|^  = 0 
ax  ay  ^ ' ax  ^ ^ ^ ' 
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where  w * | Is  negligible  by  definition  of  h(x,y,t). 

-h 

Ing  the  free  surface  condition,  we  have 


Substltut- 


au  ^ ^ ^ . 

ax  ay  9t 


(21) 


Equations  (19)  and  (20)  are  now  further  simplified  by 
neglecting  the  non-linear  terms.  The  reader  is  referred  to 
Appendix  2 for  a detailed  plausibility  argument  concerning  this 
step*. 

Two  final  simplifications  will  be  made  in  equations  (19) 

and  (20),  The  Coriolis  parameter  2Qsin(~)  will  be  linearized 

R 

by  writing  2 ..'2  sin(-^)  ~ py  where  p = 2Q/R, 

In  addition,  if  the  velocities  are  found  in  some  manner, 
then  the  free  surface  shape  can  be  obtained  by  integrating  the 
equations  (19)  and  (20)  (neglecting  the  Integrals  of  the  non- 
linear terms)  with  respect  to  x and  y respectively.  This  yields 

(sDPl  + i^  g-1^)  =X 


Where  X denotes  a known  function.  The  solution  of  this  quad- 
ratic equation  in  q is 

- D + D Jl  + — i JL 
y gP 

n = 

Ap 


* It  must  be  emphasized  that  the  argument  presented  in  Appendix 
2 is  one  of  plausibility  and  not  one  of  justification.  In 
view  of  the  desirability  of  obtaining  an  analytic  solution, 
we  neglect  the  non-linear  terms  in  the  hope  that  the  results 
will  agree  qualitatively  with  observation  and  v/ill  so  furnish 
a mathematical  description  of  the  ocean  circulation. 
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Hence 


1 + -h  XX  31+fXXX 
Ap  ])2  i\p  gp  p2 

2p 

4:  -^  < 1. 

Ap  gp  D 


nz-^x 


and  the  pressure  term  can  be  approximated  by 


- gD 


provided  the  above  inequality  holds.  It  will  be  shown  in  Sec- 
tion 5 that  the  values  of  the  constants  v/hich  are  appropriate 
to  our  problem  satisfy  this  condition. 

Hence,  the  final  equations  take  the  form 

^ « ByV  = - gD  |H£.  + MU  + (22; 

at  ax  X 

M + pyU  = - gD  ^ + MV  + (23: 

iS  + ^ = - l2£  . (2k' 

ax  ay  at  ^ 

The  boundary  conditions  are  U = V = 0 on  a land-water 
boundary.  The  wind-stress  is  prescribed  to  be 


T^x  = “ + P’  sin  wt)cos  ny,  Ty.  = 0 


where  W*,  P'  represent  the  magnitude  of  the  mean  wind-stress 
and  the  amplitude  of  the  time  variation  of  the 
wind -stress,  respectively, 

u)  is  the  frequency  of  the  wind  variation, 

n is  the  wave  number  associated  with  the  wind  dis- 

tribution. 
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One  can  consider  the  above  form  for  the  wind  as  a typi- 
cal term  In  a Fourier  series  for  a more  general  wind  distribu- 
tion. The  numerical  results  in  this  report  are  based  on  a value 
of  CO  corresponding  to  a period  of  one  year  and  n is  set  equal 
to  2it/s  where  s is  the  north-south  length  of  the  ocean. 

The  problem  defined  by  equations  (22),  (23),  (2>+)  to- 
gether with  the  boundary  conditions  and  the  wind-stress  term 
will  be  referred  to  as  the  one-layer  problem  or  Problem  1; 

("one  layer"  because  the  integration  over  z is  carried  out  over 
the  entire  depth). 

For  the  second  problem  in  which  the  density  stratifica- 
tion is  specified  as  two  constant  density  layers,  we  have  equa- 
tions (5)  - (10),  Each  equation  will  be  integrated  over  the 
vertical  coordinate,  z,  with  (5)  - (7)  integrated  over  the  top 
layer,  i.e.  , from  z = D2  + ^2  'to  z ^1  + ’ll?  (8)  - (10) 
integrated  over  the  lower  layer,  i.e.  , from  z = 0 to  z = Dg  + qg. 

As  in  problem  1,  the  non-linear  terms,  u^l^  9r)/9t  etc., 
resulting  from  the  Interchange  of  differentiation  and  integra- 
tion,are  neglected.  The  viscous  terms  are  Integrated  in  the 
same  manner  and  the  Coriolis  parameter  is  again  linearized. 

Then  the  Integrated  forms  of  (5)  - (10)  are 

— - - PyVi  + g(Di  - D2  + = -^^1  '^Ix  ■'^2x 

^ + PyU^  + g(Di  - Dg  + r,^  - Hg)  + '^ly  " '^2y  ^^6) 

au.  avi 


(27) 
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au, 

c 

at 


av 


Q ** 

- PyVg  + g(D2  + + Pini>AAU2'^'^2x"'^, 


ox 


— + PyU2  + g(D2  + r]^)-l^[hp2^2  + P ini3=MV2y"'^( 


at 


where 


®°2  . ='^2  - _ 8 


2y  oy 


Ui  = 


3x 

|)Di+tii 


9y 


■ " H ^P2^2^ 

nDi+Tli 


(28) 


(29) 


(30) 


p^u^dz, 


Vi  = 


D2+^2 

>D2+^2 


PlV^dz, 


^2  = 


P2V^^ 


V2  = 1 


UD2+T)2 

pE>2+^2 


0 


J 


P v^dz, 
2 2 ’ 


0 


'"lx’ 

'^ly 

'^2x» 

""2y 

^ox’ 

T 

oy 

T , are  the  x and  y components,  respectively  of  the  wind- 
-^y  stress  on  the  free  surface 

are  the  x and  y components,  respectively, of  the  shear 
stress  between  the  lower  layer  and  the  upper  layer  at 
the  interface, 

are  the  x and  y components,  respectively,  of  the  shear 
stress  between  v/ater  in  the  lower  layer  and  the  ocean 
bottom* 

We  specify  to  take  the  same  form  as  in  Problem  1* 
The  remaining  shear  stress  terms  are  assumed  to  be  negligible. 
The  boundary  conditions  are  U-]_  = = U2  = Vg  = 0 on  a land- 

water  boundary,  l.e.,  vanishing  mass  transport  in  each  layer. 
These  conditions  are  much  more  restrictive  than  the  boundary 
conditions  of  the  one-layer  problem  since  there  can  be  no  verti- 
cal Interchange  of  transport  across  the  interface  at  the  bound- 
aries. 

Equations  (25)  - (30),  together  with  the  boundary  condi- 
tions and  the  wind-stress,  constitute  Problem  2,  or  the  two- 
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layer  problem  (the  vertical  integration  being  carried  out  in 
two  steps). 

It  may  seem  to  the  reader  at  this  point  that,  since  we 
have  Integrated  the  equations  of  motion  over  the  vertical  coor- 
dinate z in  both  problems,  there  is  nothing  to  be  gained  by 
considering  Problem  2 in  which  the  density  distribution  is  more 
specialized  that  that  of  Problem  1.  Because  of  the  Importance 
of  this  point,  we  shall  discuss  the  significance  of  the  two 
problems  in  more  detail. 

Needless  to  say,  the  problem  of  greatest  interest  in- 
cludes the  more  general  density  distribution  of  Problem  1,  the 
four  Independent  coordinates  x,y,z,t,  and  the  full  non-linear 
equations.  The  wind-stress  components  appear  as  the  values  of 
the  vertical  shear  at  the  free  surface  z = T](x,y,t),  The  solu- 
tion of  this  problem  would,  of  course,  include  complete  inform- 
ation concerning  the  dependence  of  the  motion  on  z.  Being 
unable  to  attack  this  problem,  we  are  forced  to  integrate  the 
equations  over  z and  to  content  ourselves  with  a solution  for 
the  transport  components. 

At  first  this  Integration  over  the  vertical  coordinate, 
z,  appears  to  have  only  one  shortcoming,  viz. , a loss  of  inform- 
ation concerning  the  vertical  distribution  of  velocity.  We 
cannot,  however,  completely  afford  such  a loss  of  infornatlon 
in  the  formulation  of  the  "transport"  problem  and  some  recourse 
to  field  evidence  is  necessary.  Unfortunately,  however,  accur- 
ate observational  data  are  extremely  difficult  to  obtain.  In 
particular,  it  is  generally  held  that  the  motion  in  the  deep 
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layers  of  the  oceans  Is  negligible,  but  no  definite  conclusions 
have  been  established  to  this  effect.  It  is  because  of  this 
uncertainty  that  we  consider  the  two  separate  problems,  1 and 
2,  If  the  motion  of  deep  water  is  really  negligible,  the  pres- 
sure gradient  in  deep  water  is  also  negligible  and  the  assump- 
tions of  Problem  1 are  justified  with  the  result  that  the  thermo- 
cllne  responds  instantaneously  to  a change  in  the  free  surface 
height  provided  the  lydrostatic  pressure  assumption  is  also 
valid.  Consequently,  the  only  motion  existing  in  the  layer 
below  the  bottom  of  the  therraocllne  is  that  due  to  the  shear 
force  exerted  by  the  water  at  the  depth  z = T - d onto  the  water 
below  it.  Vertical  shear  will  extend  the  motion  to  lower  depths 
but  the  velocities  will  decay  exponentially  in  the  vertical 
direction  [l]  until  they  become  negligible. 

If  the  motion  of  deep  water  is  not  negligible,  then  we 
must  consider  Problem  2 where  no  such  assumption  is  made.  In 
that  case,  the  thermocllne  does  not  necessarily  respond  imme- 
diately to  a change  in  the  free  surface  and,  consequently,  a 
pressure  gradient  may  result.  Since  the  fluid  in  the  bottom 
layer  is  homogeneous  and  since  the  wave  length  of  the  thermo- 
cline  is  large  compared  to  the  depth  of  the  lower  layer,  a 
velocity  with  uniform  vertical  profile  is  set  up,  (hydrostatic 
pressure  being  again  assumed).  The  shear  stress,  exerted 

by  the  water  of  the  upper  layer  onto  the  surface  of  the  lower 
layer  also  causes  a velocity  in  the  lower  layer.  This  velocity 
is  not  uniform  vertically.  The  problem  including  the  effect  of 
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^2^  and,  in  addition,  the  stress  of  the  ocean  bottom  on  the 
lower  layer,  is  so  complex  that  an  analytic  solution  is  out  of 
the  question.  V/e  therefore  assume  that  the  effects  of  these 
shear  stresses  on  the  velocity  in  the  lower  layer  are  negligible 
when  compared  to  the  velocity  resulting  from  the  variation  of 
the  thermocline. 

If  the  two  problems  were  now  solved  and  the  results 
compared  with  available  observational  data,  it  might  be  possible 
to  determine  whether  or  not  sensible  deep-water  motion  exists* 

As  we  shall  see  in  Sec.  however.  Problem  2 cannot  be  solved 
by  the  methods  used  in  the  present  paper,  and  numerical  methods 
of  solution  may  have  to  be  employed. 

4.  Solution  to  Problem  1.  The  solution  to  Problem  1 
will  be  carried  out  by  means  of  a boundary  layer  technique.  For 
the  convenience  of  the  reader  who  is  not  familiar  with  this 
technique  and  who  wishes  to  follow  the  details  of  the  present 
section,  a discussion  of  boundary  layer  analysis  is  presented 
in  Appendix 

The  solution  of  differential  equations  by  boundary  layer 
analysis  can  be  carried  out  most  conveniently  if  the  equations 
are  first  put  into  non-dimensional  form.  Let  the  rectangular 
ocean  have  dimensions 

0 < X < r^,  0 < y < s (Fig.  1). 

Choose  as  a reference  length  the  north-south  dimension, 
s,  and  define  dimensionless  coordinates  x’,  y*  by 
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y = sy',  X = sx^ 


Then  the  east-west  and  north-south  dimensions  of  the  ocean  in 
non-dimensional  coordinates  will  be 

0<x*<— i = r,  0<y*<l, 

We  shall  assume  that  the  ocean  is  bounded  by  land  on 
X*  = 0,r  and  by  water  on  y*  = 0,1, 

Now  differentiate  equation  (3, 23)  with  respect  to  x and 
equation  (3*22)  with  respect  to  y and  subtract.  Substituting 
for  the  prescribed  wind-stress,  we  then  have 


0t^8x  “ 8y 


) ^ py(|S 


* Ij)  + PV 


= AA(M  . 2U) 
8X  . ay' 


- [nW»  + hr*  sin  ut]sin  ny. 


(1) 


Introducing 


and  defining 

nW»  = W,  nr'  = r, 
equation  (1)  becomes 


if  - M-)  + Bv'(^  + + flV 

s 8t''8x'  ^ av*-'  ^ P'' 


9y 


8x 


ay' 


+ 


8x»8y 


2 


a^u  _ 8^  ^ 

8x'^8y  0 y*^ 


- W [1  + a sin  T]sin  nsy' 


or 


(2) 


*i 
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jj_  _9_  p ay 
sV/  ^ ax  • 

au  -1 
" ay'-* 

w '‘ax* 

+ £ V 

w 

1 — 1 

a^v  + 

a^v 

a^u 

3- 

a'^ 

r 

* 

ax»3 

ax*  dy*^ 

<0 

CM 

<0 

ay 

- [1  + a sin  t]  sin  nsy'.  (3) 

Now,  since  the  term  (1  + a sin  T)sin  nsy*  is  of  order 
unity  } and  since  this  term  represents  the  wind  \\/hich  generates 
the  velocities,  it  is  appropriate  to  choose  a dimensionless 
velocity  which  will  also  be  of  order  unity.  Hence  we  select  a 
non-i-dimensional  term  containing  the  velocity  which  is  presumably 
of  order  one.  The  term  suggested  by  an  inspection  of  (3)  is 
pv/W  and  we  therefore  put 

V = M and  U = IP  , 

w w 

V/e  shall  drop  the  primes  from  the  x'  and  y'  coordinates 
and  work  in  the  non-dimensional  system  henceforth.  With  the 
definitions,  e = A/ps^  and  6 = w/p«  equation  (3)  becomes 

6[V^  - + Vy]  + V = e[V^^  + V^yy  - - Uyyy  ] 


- (1  + a sin  T)sin  nsy  (^-) 
where  = av/3x,  (V^  - Uy)^  £ - a^U/ayde  , etc. 

If  we  non-dimensionalize  the  momentum  equations  (3,22) 
and  (3.23)  and  the  continuity  equation  (3.24)  by  means  of  the 
above  definitions,  we  must  introduce  a new  parameter  9 and  a 
variable  H defined  by 


9 = 1!^,  H = 


w 


* As  will  be  seen  later,  we  shall  choose  a to  be  0.2, 
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The  equations  become 


ns6  ^ - nsy  V + ©iS  = nseAU  - (1+a  sin  t)cos  nsy  (5) 

dx  3x 


nsb  ^ + nsy  U + 6 = ns  e AV 

d'z  9y 


and  (3i8)  becomes 


M + M = ^ M , 
ax  ay  “ ° 9t  * 


Attempts  to  solve  equations  (5)  to  (8)  in  closed  form 
were  unsuccessful.  We  therefore  resorted  to  seeking  solutions 
by  a perturbation  expansion  in  the  parameter  6. 


U = U + 6Ut  + + *,, 


V = Vq  + + ... 


H = Hq  + + 6 H2  + ...  . 

Our  formal  procedure  is  to  regard  the  coefficients  U^,  U,  etc,  j 
* 

as  coefficients  in  a power  series  in  b. 

Let  us  substitute  the  expansions  into  equations  (4), 
(5),  (6)  and  (7),  We  have 


6[Vn,  + 6V,.,  + ...  - - 6D 


ly  ■ ••• 


+ y [Uqx  + ‘Ulx  + •••  ^nv  + + ...] 


oy  ^ «’ly 


+ Vq  + 6V;j_  + ...  - + ... 


^oxyy  ^^Ixyy  “ ^oxxy  ~ ^^Ixxy  “ • • • 


Igyyy  » ^^lyyy  (l  + 0 sin  T)sin  nsy  (8) 
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au.  ^ ^ 

ns6[ — 2 + 6 — i + ...  ] - nsyCVp,  + &Vn  + « . , ] 

an. 

+ ©[■g^  + a + • ♦ • ] - nseA  [Uq  + 67^  + • • « ] 


- (1  + a sim;  )cos  nsy 


(9) 


av 


f!c 

ax 


+ 6 


3V, 

+ 6 

-J=  + 

3t7 

...  ] - 

(-  nsy[UQ 

+ 6 U-j^  + ...  1 

+ ©c 

3H„ 

3Ht 

0 ^ 

9y 

b 

ay 

+ • . . ] 

= nseA  [ Vq  + 6V^  + , 

...] 

3U^ 

+ ... 

av 

aVi 

+ b 

r 9H^ 

+ ...  = - 6L — 2 + b 

3Hi 

3x 

9y 

9y 

d'v 

3t 

(10) 


(11) 


If  we  regroup  each  of  these  equations  so  as  to  combine 
the  coefficients  of  each  power  of  6 y v/e  have,  upon  retaining 
terms  in  6°  and  6 only: 

v[U  + V 1 + V - efV  + V - U - U 1 
ox  oy-*  0 '"'-''oxxx  ^ '^oxyy  ^oxxy  ^oyyy* 

+ (1  + a sin  -c)sln  nsy  j'  + I - Uoy,,  + y[Ujjj  + Vj^y] 

* h - - hxyy  - ''ixxy  ' °lyyy^  } =0  C12) 

0H  1 

-nsyVQ  + © - nseAUQ  + (1  + a sin  t:)cos  nsy  p 


+ i ns 


3U 

I 


an. 


■j  nsyU^  + © - ns  eAV 


9Hn  1 

2 - nsyVx  + © - nseAU^  ^6  + ...  = 0 (I3) 

1.  J 3Hn 

„ >+*<ns  — 2 + nsUT  + © — — i-nseAVi^-tjM  = 

J L 3-^  ^ T,,n^ 


0 


£Ho  ^ilo 

3x  3y 


j 

fa^  9^  mJ 

M 

[_  8x  3y  3t  j 

J ilk) 
(15) 


Setting  each  of  the  coefficients  of  6 equal  to  zero  we 
have  as  the  zero  order  equations  for  (12)  and  (l5) 


All- 101 


ko 


^^oxxx  ^oxyy  “ ^oxxy  " ^oyyy  ° T^)sin  nsy  (16) 


UOX  = 0 


(17) 


and  the  boundary  conditions 


Uq  = Vq  = 0 on  X = 0,  X = r. 


(18) 


V/ith  the  particular  wind  distribution  prescribed  we  will 
also  be  able  to  satisfy  the  additional  boundary  conditions 

Vq  = = 0 on  y = 0,1.  (l8.a) 

We  shall  proceed  to  solve  equations  (l6),  (17)  together 
with  the  boundary  conditions  (18),  (l8. a)  for  the  velocities 
Uo  and  Vq. 


Define  a stream  function 


V = 

0 8x^ 


TI  - 


(19) 


SO  that  (17)  is  satisfied  identically.  Then  (l6)  can  be  written 


eM^|)  = (1  + a sin  T)sln  nsy 

where  AA(  ) is  the  blharmonic  operator  + 2 ^ + 

L.  . 3x^  3x^ay2 


(20) 


Equation  (30)  is  similar  to  the  one  solved  by  Munk 
and  I'^nk  and  Carrier  [6].  In  the  present  case,  however,  the 
non-dimensional  time,  t , appears  as  a parameter,  so  that  our 
problem  corresponds  to  a quasi-steady  problem. 

Equation  (20)  together  with  the  boundary  conditions 
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= 'vj)  =:0  onx  = 0,r 


^ = 4>yy  =0  on  y = 0,1 


(20. a) 


can  be  solved  for  by  applying  the  boundary- layer  technique 
to  the  boundaries  x = 0,r.  The  solution  is 

{“1/3 

-X  + r - 


l/r,  xV3e 

+ [ ( e -r)cos( 


-1/3 


•)  + 


1/3 


X Vie 


J./3 


X e 


-l/3> 


+ ( '/S  e sin( ) ] e ^ 

\/3 

From  (19)  and  Vq  are  found  to  be 

r 

Uq  = - ns(l  +a  sin  'r)cos  nsy-j-x  + r - + e^/3, 


(x-r)e 


(21) 

-1/3 


+ [(  e^^^-r)cos(i)&|I )+(\/3e  ^ L)sin(ii^ )]  e 

^ l/3  ^ 


-1/3^ 


(22) 


Vq  = (1  + a sin  x)sin  nsy 


, ,2re 

+ [C0S( 7^ ) +( 


{■> 

-1/3 


(x-r)e‘^^3 


+ e 


V3 


-1/3  - 11. 

-)]e 


-1/3 


-V3)sin(-'^^  2 


I 


(23) 


The  zero-order  equations  derived  from  (I3)  and  (l4)  are 


* The  problem  defined  by  equations  (20),  (20. a)  is  solved  in 
detail  in  Appendix  5 by  means  of  the  boundary  layer  technique. 
The  method  used  in  the  remainder  of  this  paper  is  described 
in  detail  in  that  section.  Munk  and  Carrier  [6 ] used  this 
method  for  solving  the  steady  problem  in  a triangular  ocean. 
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= nsyVo  + useAUq  - (1  + a slnT)cos  nsy  (24) 


®^oy  " “ + nssAV^. 


(25) 


Solving  for  Hq,  we  have,  (neglecting  terms  of  order  e) , 

1/3 

QHq  = (1  + a sinT)(cos  nsy  + nsy  sin  nsy)(  - x + r - e ) 


+ (1  + a sin  T)nsy  sin  nsy  ■<( 

+ [(e^'^3  . r)oos(iLVl£^)  + 


1/3  (x-r)e 
e e 


-1/3 


+ ( \/3 


- -^)sin(iL^i 


-1/3 


/3 


Ole 


-1/3] 

xe 

->  * 


J 


(26) 


First-Order  Solution 

From  equations  (12)  and  (l5)  the  terms  of  first  order 
in  6 are  found  to  be 


'"ixyy  ' “ixxy  ’ "'lyyy  ’ ' ''l  = t'^ox  ’ "oy  ’ U7) 


°lx  " ’'ly  = - Ho^ 


(28) 


The  boundary  conditions  are  again  = 0 on  x = 0,r, 

In  (27)  and  (28)  the  right  sides  of  the  equations  pro- 
vide the  driving  term  as  did  (1  + a sin  T)sin  nsy  in  the  zero- 
order  equation.  We  shall  proceed  with  the  solution  by  means  of 
the  boundary  layer  technique* 

For  the  interior  solution  we  assume  that  the  functions 
are  smooth  and  hence  that  the  derivatives  are  of  the  same  order 
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of  magnitude  as  the  functions  themselves.  The  terms  multiplied 
by  e may  therefore  be  neglected. 

Let  us  rewrite  equations  (22),  (23),  and  (26)  as  the 
sum  of  two  parts  - one  part,  with  subscript  i,  having  the  same 
order  of  magnitude  throughout  the  domain  (the  "interior  solution"); 
the  second  part,  with  subscript  b,  sensibly  large  near  the  bound- 
ary and  negligibly  small  in  the  interior,  (the  "boundary  layer 
contribution" ) 


1/3 

= - ns(l  + a sin  t:  )cos  nsy(-  x + r - e ) 


= - ns(l  +a  sin  t)cos  nsy 


,1/3  ,(x-r)  -1/3  ^ 


4.  [ _L)sin(,^4f.,,,)]e  2 

2 1/32 


-ly^' 


Ob 


V.  = -(l+a  sin  T)sln  nsy 

01  -l/"^ 

f (x-r)e"^^^+  [cos(2Li5i ) + 

(1  + a sin  T)sin  nsy”3e  2 


o -1/3 

+ (2re \/3)sin( 

/3 


\/3 


X 


-1/3  - xi 


-1/3  n 


)]e 


I/O 

= (1  + a sin  t)(cos  nsy+nsy  sin  nsy)(-  x + r - e '^) 

[ 1/^  (r-r)e”^^^ 

©H^b  = (1  + a sin  T)nsy  sin  nsyi  e '-^e  + 


+ [(6^^3.^)^Pg(j^/3.f.^^^)^(|/3el/3,  _x,)sin(iL^i-.^^]e  ^ 

2 Vi  2 


-131 


) 
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We  expect  the  boundary  layer  thickness  to  have  the  same 

order  of  magnitude  in  the  higher  order  solutions  as  in  the  zero- 

1/0 

order  solution,  viz»  , e ' Thus,  in  order  to  find  the  first- 
order  interior  solution,  we  neglect  all  the  terms  with  subscript 
b since  they  are  negligible  in  the  interior.  Thus  immediately, 
the  interior  portion  of  is  known  and  is  (from  (27)) 

"'ll  = - f''olx  - “oly  - yHoiW 

= ° (-X  + r - e^^^)[cos  nsy  + (nsy-On^s^)sln  nsy] 

(29) 

From  (28)  and  (29)  the  interior  portion  of  U^,  can 

be  computed  directly,  giving 

= - .jS .££.g.  J.  [ - ^ + x(r  - e^'^^)]  [2nsy  sin  nsy  + 

2 2 2 

+ (n  s y + ©n-^s-^  + 2)cos  nsy  ] + C^(y,T) 

where  C^(y,T)  is  arbitrary  and  must  be  evaluated  by  applying 
the  boundary  conditions  to  the  complete  solution,  l.e. , inter- 
ior solution  plus  boundary  layer  contribution. 

Before  proceeding  with  the  boundary  layer  analysis  we 
can  simplify  equation  (27)  to  some  extent.  Near  x = r. 

Vox  = 0(e"^/3)^  ^ 0(e^/3)^  and  = ©“^0(e^/3),  N^ar 

X = 0,  = 0(e“^^3)j  Uay  0(1),  and  Hay  = ©“^0(1),  Thus  in 

each  case  we  are  justified  in  using  only  the  contribution  of  the 

-2/3  -1 

V term  provided  e >>1  and  e ■^'-^>>  © , As  will  be  shown 

later,  when  the  appropriate  dimensional  constants  are  substi- 
tuted, the  error  involved  in  neglecting  the  other  terms  is 
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^5 


extremely  small.  Thus  for  all  practical  purposes,  equation 
(27),  near  the  boundaries  can  be  written 


Ibxxx  ^Ibxyy  “ ^Ibxxy  ” ^Ibyyy  ^ ^oxt 

J -I/3  (x-r)e 

= ft  cos  T sin  nsv  1 e “*6  + 


= a cos  T sin  nsy  1 e -^e 

— 1/*^ 

-2/3  -1/3  . xe"^^^ 

sin(iL^|l_)]e  2 

1/3 


(31) 


Near  x = 0,  the  inhomogeneous  contribution  which  contains  the 


term  e 


(x-r)e 


-1/3 


can  be  neglected  since  its  effect  is  felt  only 


near  the  eastern  boundary,  i. e. , near^^^=  r.  Similarly,  near 
X = r,  the  terms  multiplied  by  e“  can  be  neglected.  Thus 

for  the  region  near  x = 0, 

^ ^^Ibxxx  ^Ibxyy  “ ^Ibxxy  " ^Ibyyy  ^ " ^ib 

= a cos  T sin  nsy[(re  ^^^-2e'’^^^)cos(.^ ) + 


-2/3 

+ ^ sinC 

V/3 


,/r  -1/3  . 2LL_ 

)]e  2 


-1/3 


(32) 


Now  suppose  the  x coordinate  is  stretched  by  substitut- 
ing X = (k  > 0),  Then  (32)  becomes 


el-3ky  4.  e^“^V  - e^**2kg  - eU  - V 

lb5^^  Ib^yy  Ib^^y  ibyyy  lb 


= a cos  'V  sin 


nsy  [(re~^^^-2e~^^3)cos(.5~ 


k-1/3 


V3  2 


) 
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The  term  of  highest  order  derivative  in  ^ is  matched 
with  the  remaining  largest  term  in  the  equation.  Hence,  we 
formally  match  with  Then  k = I/3*  and  the 


equation  becomes 


^Ib^^^  ” ^Ib  “ cc  CIOS  T sin  nsy[ (re*’^^^-2e’*^'^^)cos(.^Ji5)  + 


r£!2  sln(ii5)] 
VI  2 


(33) 


The  term  can  now  be  expanded  in  an  asymptotic  series 
in  e and  only  the  first  terms  will  be  kept.  Since  the  inhomo- 
geneous term  of  (33)  contains  only  exponential  and  trigonometric 
functions,  let  us  try  a solution  of  the  form 

= a cos  V sin  nsye“^^^/v°  cos(-i^)+V°  sin(-Lk!3) 

L ^ J (3^) 

where  and  are  the  first  terms  of  asymptotic  expansions 
and  are  to  be  determined. 

If  as  given  by  (34)  be  substituted  into  (33)  and 
s iri  X \/o 

if  coefficients  of  ^qs  (— ^^)  be  equated,  two  simultaneous  differ- 
ential equations  with  constant  coefficients  result* 


- 1 v°  - 3 v°  ^ + v° - 3JZ3  + ly^. 

2 1^  2 l^C  UU  2 U 2 


3 yO  3 yO  _ 3 yO  . 3 yO  . yO  _ j. 

~ — •''la  2 1^  2 ^IKK  ^ 


V°  = r - 2e^^^ 
Vi^^  - r - 


(35) 


(36) 


* The  fact  that  k = I/3  indicates  that  the  thickness  of  the 
boundary  layer  is  of  the  same  order  of  magnitude  in  the  zero 
and  first  order  solution,  as  was  anticipated. 


^7 


f ! 

^ i 

f 


Particular  solutions  of  (35),  (36)  are 

V3  „ ,1/3  ^ 

V°=! v°  = 3-e  . -.r 

3 ^3/3 

The  homogeneous  solutions  may  be  derived  by  letting 


r°  = , V°  = Be^^  , 


Then  (35),  (36)  become 


I X - I + x^)  - B [i^  X-  X^  ] = 0 


(37) 


A(3^  X - X^)  + B [-  |x  - |x^  + \^]  = 0.  (38) 


Hence,  since  the  determinant  of  these  two  simultaneous 


equations  must  vanish,  we  have 


(>.5  - - I X)®  (X-  x^)^  = 0. 


(39) 


The  roots  are 


x = 0,0,  i-1^  , , )/31,  - v/stClfO) 


Then, 

1/3  3-^Vll  v/315  3 -/3I,  - 

v°  = J =1  ^ + A^e  2 Age  -fA  e 2 A^e  + A^ 


0 liiSl,  \/315  >^5  .y5U 

7?  = ^ ^ + B^e  2 ^ + Bpe  +B,e  2 ^ +b.  e % 

1 3 \/3  5 ^ b -+ 

Hence,  from  (3^) 


^Ib  ~ e "2/3  ^ A^e^^^+  A^^e^'^^^Aj] 

^ sin  (-^^)  [ 3^'"^-.  ^ * B^e 

(>+l) 


where  we  have  set  = B ^ = 0 since  the  contributions 

of  the  terms  with  those  coefficients  do  not  tend  to  zero  as 
^ > 00  , 

When  (37)  and  (38)  are  used  to  get  a relationship  be- 
tween the  and  the  B^,  then  the  final  form  for  near  x = 0 


is  found  to  be 


V = a cos  T sin  nsy  e 
lb 


-2/3 


1/3 


02)cos(i^) 


1/3 


- r 


3 V/3 


K + C^)  sin(V^)|e"^'^^  (42) 


where  and  are  arbitrary  functions  of  y and  1 and  must  be 
found  by  applying  the  boundary  conditions  to  the  complete  solu- 
tion. 

In  a similar  manner,  if  we  make  the  following  two  sub- 
stitutions for  the  right  (eastern)  boundary 

hn 


(x  - r)  = e 


Vib  = ot  cos  t;  sin  nsy  e 


-l/3„ne 


h-1/3 


e^"  [VO+...I, 


we  find  that  h = 1/3  and 


"I/3 

V.,  = a cos  T sin  nsy  e [ U + An  (y , m)]  e ^ , 
lb  3 

We  have  used  the  fact  that  — > 0 as  t]  — > - co . (As  stated  in 

the  appendix,  rj  -->-<"00  v;hen  the  boundary  on  the  right  is  under 
consideration,  since  the  boundary  layer  solutions  must  become 


* The  same  remark  applies  to  the  value  of  h as  previously  made 
for  the  value  of  k. 


1,-TOl 


h9 


negligibly  small  as  the  distance  from  the  boundary  increases, 
i,  e.  , as  T)  or  X decreases* 

If  the  three  contributions  (29),  (42),  (43)  to  the  com- 
plete solution  for  Vj_  be  added,  the  final  form  for  is 

„ a cos  Tx  r/  2 ^ 2 2n  . T 

= (-X  + r - e ) [(y  ns  + 9n  s )sin  nsy  + y cos  nsyJ 

-L  /-s 


+ a cos  T sin  nsy  e 

-1/3 


{(1-.-  X + C2(y))cos( 

L 3 . n 


/3 


-1/3 


-1/3  1 - xe~^'^3 


+ X + G^(y))sin(---i:-|^ ) je 

. • . r -1/3  , , 

+ a cos  T sin  nsy  e e + A^(y)  re 


2 

-1/3 

(44) 


U, 


By  means  of  the  continuity  equation  we  then  find 
^ [2nsy  sin  nsy  + (y^n^s^  + 2 + On^S'^)cos  nsy] 


[-  ^ + x(r  - 6^^^)  ] + C^(y)  - ^ 

TA  1 (x-r )e”^'^^  T (x-r)e“^^^ 
- a cos  T ns  cos  nsy  [At-  A + ] e 

■^3  3 


^insy  sin  nsys^/Set^-':')'!'^^ 


8A,  (x-r)e”^'^^ 

- a cos  T sin  nsy  — — i e - &,-9J5.§.X 

9y 


•[(re^'^^  - 2e2/^)cos( 


9 

i/=7  “1/3  1/^  • y-  -1/3  - X e 

)]e  ■ 2 


nsy  sin  nsy  • 


-1/3 


- a cos  T ns  cos  nsye 
2/3  . 1/3 


/3 

-2/3j-(.2_xe^  _ 


+ ^)cos(ii^ 


-1/3 


•)  + 


+ (- 


„ i/d  \/^  -1/3  - X g- 

- -12^.  - 2^_)sin(2LVl£.„)]0  2 


-1/3 


v/3  3 3 

e-^3  , f 


+ a cos 't 


— -fel  (C2  ^3)  sin  nsy  cos(— 

^ I n - 


/3 


,1/3 


"1/3  "1  - xe' 

+ (C^  - \/3  C2)sin  nsy  sin  (i-O^-)  > 


■1/- 


■)  + 
(45) 


1 
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The  arbitrary  functions  of  y can  be  evaluated  by  means 
of  the  boundary  conditions  = 0 on  x = 0,r.  V/e  have 

2/3  222 

sin  nsy  C2  = [(nsy  + 6n  s )sln  nsy  + y cos  nsy]  (4-6) 

-2/3  222 

sin  nsy  A-j^  = -- — [(nsy  + ©n  s )sin  nsy  + y cos  nsy]  (47) 

= -S  .Q.P.§ .S.>|^[2nsy  sin  nsy  + (y^n^s^  + 2 + 9n3s3)  • 


* cos  nsy][^  - re^^^  + + l)cos  nsy  “* 


(48) 


sin  nsy  C = [ 2e^/^(y^ns  + ^ + ©n2s2)(re^/3  . - e2/3)  + 

^ \/3  © L ns 

+ ~r  — ]sln  nsy  + (5y  cos  nsy  — ^ sin  nsy)re^'^^  - 

3 3 ^ ^ ns 

2/3  2/3  o I 

- (9y  cos  nsy  - -iS  sin  nsy)e  '^+(re  ' -e) (y^ns+©n^s^)sin  nsyk 

ns  J 


The  first-order  contribution  to  H can  be  found  from 
equations  (12),  (13)»  The  first  order  equations  are 

9U  3H, 

n e O o’Tr  V J.  Q X _ v,t,^ATT 


- nsy  V.  + © — i = nseAU- 
1 3x  1 

9H 

+ nsy  % + © ^ = nssAV, 
oy  1 


from  which  is  found  to  be 


h-- 


nsa  cos  t 


•|[(0ns  + y^)cos  nsy  + (y^ns  + y©n^s^)sin  nsy] 


• [kx^  + r2)  + (e^/3  _ ^ gl/3)  ] + _©, 

^ ns 

g2/3(©^  ^ ^^^y  sin  nsy  ^ cos  nsy^\^ 

3 ns  n^s2  J 


cos  nsy  + 
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+ nsy  sin  nsy  -) 


3 1/3  3V^- 


+ [-(r  _ gi^ 
'■3  3 


,2/3 


)..1^3eos(-V-3s“^/^  sI/3 


- M. 


2 

-1/3 


) - -^i(C^+  \/3C^)cos( 


x\f3i 


+ (C^  - 1/3  C2)sin(  )]  pe  2 + g„-|QS  ..T.  nsy  sin  nsy 


I (X  - r - e^/3,  -1/3  ^ |.^(x-r)e 


-1/3 


(50) 


The  terms  U-,  and  do  not  satisfy  the  boundary  con- 

au.  ^ 

ditions  = 0 on  y = 0,1*  We  must  recall  that  these 

boundary  conditions  were  chosen  rather  arbitrarily  as  being 
plausible  ones  for  the  type  of  wind  distribution  specified,  and 
the  y dependence  of  the  zero-order  solution  was  accordingly 
chosen  as  sin  nsy«  VJe  cannot  expect  such  a y dependence  to 
satisfy  all  the  conditions  for  each  set  of  equations.  The  fact 
that  U-j^  and  V-j^  do  not  satisfy  the  boundary  conditions  does  not 
seem  to  be  very  serious  since  we  do  not  really  know  what  con- 
ditions are  appropriate. 

If  we  next  consider  the  equations  resulting  from  equat- 

2 

Ing  the  coefficients  of  6 to  zero,  we  obtain  from  (8)  and  (11) 


^^^2xxx  ^2xyy  ~ ^2xxy  " ^2yyy  ^ " ^2  “ " ^ly"y^l^T; 

^2x  + ^2y  = - ^It  • 

In  the  boundary  layer,  near  x = 0,  is  of  order  e“^ 
Thus  we  can  expect  V2  to  be  of  order  e"^  in  that  region.  By  a 
similar  argument,  we  can  expect  to  be  of  order  e“^/3^  to 
be  of  order  e ^'-’,etc.  If  we  therefore  write  out  the  series 
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V = Vq  + 6^2^  + 5 V2  + ^^^3  ^ ^ ^If  ■*■••* 
we  have  In  terms  of  orders  of  magnitude  near  x = 0, 

V = 0(e“^'^^)  + Q^g-l/S)  + 

+ 0(e“^'^^)  + fc^e-l  0(e“^^)  + , 

-1/3 

or  factoring  out  the  0(e  ) , we  have 

V = 0(e“^^^)[l  + (5e“^'  ^ 


The  perturbation  scheme  may  be  expected  to  be  valid 
-1/3 

provided  6e  ^1,  We  can  expect  a fairly  good  approximation 

from  only  the  first  two  terras  provided  the  more  stringent  con- 

-1/3  -1/3 

dition  6e  '‘^<<lis  Imposed.  If  6e  = l/5>  the  error 

involved  in  neglecting  the  third  term  is  no  larger  than  of 

the  first  term. 

-I/3 

For  yearly  variation  of  the  wind,  be  ~ 1/6,  Hence 
we  shall  keep  only  the  first  two  terms  of  the  series.  It  should 
be  noted  that  a determines  the  magnitude  of  the  effect  of  the 
perturbation  but  it  has  no  bearing  on  the  validity  of  the  ex- 
pansion. 


Numerical  Example 

In  order  to  discuss  the  above  solution,  we  shall  pre- 
scribe numerical  values  for  the  constants  of  the  problem.  Let 


r-j^  = 6.  5 X 10  cm 
g 

s = 5 X 10  cm 


P = 2 X lO^^^cm”^  sec"^ 

D = 5 X 10*^cm(C  = 200m.  , d = 600m. ) 
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_ 7 2-1  “7  -1* 

A = 5 X lO'cm  sec  ^ u = 2 x 10  sec 

t]  = 2Tt/s  V/»  =0.6^  gm  cm“^  sec"^  . 

The  magnitudes  of  r^,  s,  A,  D correspond  roughly  to  the 
Atlantic  Ocean  parameters.  The  value  of  p is  chosen  so  as  to 
give  the  best  approximation  to  the  Coriolis  parameter  in  the 
latitude  of  Cape  Hatteras.  The  equality  n = 2%/ s corresponds 
roughly  to  the  east-west  components  of  the  trades  and  the  west- 
erlies. The  value  of  w corresponds  to  yearly  frequency  of  the 
wind  variation  and  W*  = 0.65  gm  cm"  sec"'^  is  the  value  used  by 
Munk  [5]  for  the  wind  stress. 

Then  the  dimensionless  constants  have  the  values 

6 = = 2 X 10"^  ns  = 2% 

e =-A-  = 2 X 10-^  9 

ps^  p S-^ 

r = 1.3 

Also  r*  has  been  chosen  so  that 

a = 0. 2, 

The  results  for  this  numerical  example  are  shown  in 
Figs.  2-6. 

In  Fig,  2 the  non-dimensional,  north-south  component,  V, 
of  the  mass  transport  is  plotted  against  x'  near  x’  = 0 for  the 
value  y'  = 0,25.  The  region  of  large  V corresponds  to  the  Gulf 


* Corresponding  to  an  annual  period  for  the  wind  fluctuation. 
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stream  and  the  section  adjacent  to  the  Gulf  Stream,  with  nega- 
tive V,  corresponds  to  the  off-shore  counter-current* 

For  the  Gulf  Stream,  the  extreme  values  of  V are  in 
phase  with  the  extreme  values  of  the  wind.  However,  for  the 
points  between  the  maximum  and  minimum  values  of  wind  strength, 
the  transport  lags  behind  the  wind. 

During  one  cycle  of  wind  variation  the  following  result 
is  found.  The  transport  and  wind  both  have  maximum  values  at 
T = %/2,  Immediately  after  t = it/2,  the  wind  begins  to  decrease. 
The  transport  also  decreases  but  it  lags  behind  the  wind.  At 
T = It  the  wind  has  reached  its  mean  amplitude  and  the  lag  of  the 
transport  is  greatest,  viz,,  an  interval  of  9 days  elapses 
between  the  time  the  wind  reaches  its  mean  amplitude  and  the 
time  at  which  the  transport  reaches  its  mean  amplitude.  After 
-c  = 7t  , the  transport  begins  to  gain  on  the  wind  until  at 
t:  = 3n:/2,  the  two  are  again  in  phase.  The  wind  and  the  trans- 
port now  begin  to  Increase  and  the  transport  again  lags  behind 
the  wind.  The  maximum  lag  is  reached  at  t = 2ti;  at  which  point 
the  transport  begins  to  catch  up  to  the  wind.  They  are  in  phase 
again  at  t = 5^/2,  This  cycle  is  repeated  Indefinitely* 

The  discussion  presented  here  is  based  on  the  ass\imption 
that  the  first  two  terms  of  the  series  represent,  in  a sufficient-' 
ly  accurate  manner,  the  complete  solution.  One  result  of  this 
assumption  is  that  transport  reaches  its  maximum  value  at  t =%/2, 

* It  is  shown  later  that  the' Value  9 days  is  independent  of  the 
specific  value  of  the  frequency  for  slowly  varying  winds. 
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The  perturbation  contribution  vanishes  at  that  instant  since 
its  coefficient  is  cos  Thus,  no  matter  what  the  value  of  6 
(essentially,  the  frequency),  as  long  as  it  lies  within  the 
limits  necessary  for  the  validity  of  the  above  method  of  solu- 
tion, the  maximum  value  of  the  transport  will  occur  at  t = mu/2, 
m = 1,  5}  9 ...  , and  its  value  is  given  by  1 + a times  the 
steady  transport  value. 

The  interval  of  9 days  between  the  tiro.e  at  which  the 
wind  reaches  its  mean  amplitude  and  the  time  at  which  the  trans- 
port reaches  its  mean  amplitude  is  also  independent  of  the  fre- 
quency, To  show  this  let  Vq  = (1  + a sin  m)Q  and  = oL  cos 
Then  V = (1  + a sin  r )Q  + 6a  L cos  t.  Since  the  mean  value  of 
the  transport  is  V = Q,  we  can  find  the  time  at  which  this  occurs 
by  setting 

(1  + a sin  t)Q  + 6aL  cos  x = Q 


or 


tan  T 


Since  t is  small,  we  can  write  tanx  ~ x and  therefore 


X * 


Q 


Substituting  x = wt  and  6 = w/Ps,  we  have  finally 


t 


which  is  independent  of  frequency  and  a , 

It  is  apparent  from  Fig,  2 that  the  out-of-phase  effect 
is  of  relatively  greatest  importance  in  the  counter-current 
rather  than  in’  the  main  stream.  The  graph  shows  the  various 
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effects  only  up  to  the  eastern  edge  of  the  counter-current  at 
X*  = 0.1.  For  x’  >0.1  only  the  mean  position  of  the  transport 
is  plotted  since  the  deviations  from  this  mean  position  are  very 
small. 

Near  the  eastern  boundary  of  the  ocean  (Fig,  3)  and  in 
the  counter- current  region  (Fig,  2),  the  absolute  magnitude  of 
the  extreme  values  of  the  transport  (which  is  now  negative)  are 
also  in  phase  with  the  extreme  values  of  the  wind  and  the  trans- 
port lags  behind  the  wind  at  all  other  times. 

Figures  4-,  5?  and  6 show  surface  contours*  for  the 
southern  half  of  the  rectangular  ocean  for  t = 0,  n:/2,  n,  3V2. 
The  contribution  of  is  very  small  throughout  the  ocean**  and 
has  therefore  been  neglected.  Thus  the  graphs  for  t = 0 and 
X = % coincide.  This  result  is  based  on  the  assumption  that  D 
is  500  meters  in  thickness.  If  D were  increased  the  above  re- 
marks would  be  even  more  appropriate.  If  D were  decreased,  the 
contribution  of  the  perturbation  term  would  be  larger  and  we 
would  therefore  have  to  account  for  it.  The  value  of  the  first- 


♦ If  v/e  define  the  thermocline  as  the  surface  at  z = T - d/2, 
then  the  contour  lines  of  Figs,  4-,  5%  and  6,  multiplied  by 
-200  represent  the  deviation  of  the  thermocline  from  its 
equilibrium  position  at  z = - C - d/2  = - d, 

**  If  for  any  of  the  variables  the  magnitude  of  the  coefficient 
of  6 in  the  perturbation  solution  is  of  the  same  order  as  that 
of  the  zero-order  term,  the  coefficient  6 = 0,002  renders  such 
a correction  negligible.  Throughout  the  present  example,  the 
only  sizable  contribution  of  the  out-of-phase  term  is  found 
in  the  north-south  transport  V in  the  boundary  layer  where 
the  function  V Increases  by  order  e~l/3.  However,  Hq  and  H]_ 
have  the  same  order  of  magnitude  throughout  the  ocean  so  that 
the  first-order  correction  can  be  neglected  throughout. 
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order  velocities  would  also  be  altered  when  © is  changed,  V/e 
shall  consider  several  values  of  © when  we  discuss  the  deflection 
of  the  thermocline  in  the  steady  two-layer  ocean. 

The  mean  mass  transport  of  the  Gulf  Stream  (corresponding  to 
the  steady  problem)  is  26,6  x 10^  metric  tons  per  second  as  com- 
pared to  Munk's  value  [5]  of  36  x 10^  and  the  observed  value  of 
72-80  X 10  metric  tons  per  second.  Munk  [5]  used  the  east-west 
component  of  an  empirical  wind  system  and  the  discrepancy  is 
therefore  due  to  the  difference  between  the  two  wind  systems, 

At  the  time  of  maximum  (minimum)  wind  the  transport  is  20^o  higher 
(lower)  in  accord  with  the  remarks  made  previously  in  this  sec- 
tion. In  the  counter-current  the  steady  mass  tranport  is  *+,6l  x 
10  metric  tons  per  second. 

The  difference  betv/een  the  computed  and  the  observed 
values  is  not  surprising  when  one  considers  the  many  idealizing 
assumptions  made.  Such  features  as  the  straight  coast  lines, 
the  simplified  theory  of  turbulence  used,  the  neglect  of  the  non- 
linear terms,  and  a more  realistic  stress-effect  of  the  wind  on 
the  water  could  well  change  the  quantitative  results  by  a factor 
of  two  or  three. 

The  problem  as  stated  and  solved  by  the  above  method 
gives  no  sensible  east-west  variation  in  the  position  of  the 
Gulf  Stream,  but  a careful  investigation  of  the  eastern  boundary 
of  the  Gulf  Stream  shows  a very  small  narrowing  of  the  stream# 

How  well  such  a result  agrees  with  field  evidence  is  uncertain 
since  our  solution  yields  no  inshore  counter-current. 
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It  would  be  interesting  to  ascertain  how  well  our  pre- 
dicted results  agree  with  observation;  specifically,  if  the  mass 
transport  of  the  Gulf  Stream  responds  as  indicated  to  variations 
in  the  wind  and  if  the  lag  of  the  transport  is  Independent  of 
the  frequency. 


5«  Methods  of  Solution  for  Problem  2.  The  equations 
(3*  25)  “ (3»  30)  are  non-diniensionalized  below  in  order  that 
boundary  layer  theory  may  be  employed.  Using  the  arguments  of 
Section  for  the  method  of  non-dimensionalizing,  we  have 


X 

y 

1 


= sx’ 


= sy', 
= U)  t« 


\ = 


Vn  = 


= 


U,  = 


u.  = 


H.  = 


= 


If 

IT’ 

w 

Ijf 

w ^ 

Tjf 

w ^ 

"w  ' 

2 2 

P S PpTjp 

W ' 


gW 


ng(D2_  - D2) 

p ^s3 

i - ^^2 

^ * 

rs^ 

e = -4.  , 

O i 

ps-^ 

*>  = I 

P s * 


a = W » 


a = 


b = 


- Pi 


^2 

?2  - Pi 


Then  equations  (3*25)  - (3*30)  become 
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ns6  ^ - nsy^i  = - © S - [Hn  - aHp  ] + 

0T  9x  a 1 8x 


ns5 


av 

0T 


dx  a “J-  •'  0x 

+ nseAU-|_  - (1  + a sin  t:)cos  nsy  Cl) 

1 + nsyUi  = -eS-n|b[H^-  aH^  ] ^ + as^Vj_  (2) 


9Ut  avi  a , 

^ ^ = -j,4_  [H.  - aHp  ] 

9x  9y  ° a-r  2 -' 


(3) 


au, 


- yVp  = - (H3_  + bH^)  - XHp  A (H^  + bHp)  +eAU2(^) 


3t; 

av. 

0 5^  ^ yll2 


9x 


(H]_  + bH2)  - XHg  ^ (%  + bH2)  +eAV2(5) 

a Up  av2  aiip 
ax  ^ ay  ^ ax  * 


(6) 


Steady  Wind 

Let  us  first  treat  the  case  of  a steady  wind,  i.e., 
a=  0 and  9^ at  = d/d%  = 0,  and  let  us  assume  that,  in  the  case 
of  steady  motion,  there  are  no  velocities,  and  hence  no  horizon- 
tal pressure  gradient,  in  the  bottom  layer.  Equations  (4  ) - 
( 6.  ) are  then  satisfied  immediately  by 


Ug  = V2  = 0,  H2  = - i 


(7) 


and  equations  (1)  - (6)  become 


nsyV]_  = 
nsyU^  = 


a r nij  _L  nsX  j2 
a X 


[ ©H-i  + — ^ ] + nseAU,  - cos  nsy 

2b  1 


- A [ ©Hn  + ] + nseAV, 

0y  ^ 2b  1 i 


(8) 

(9) 


+ = 0 
9x  9y  * 


(10) 
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Differentiating  (8)  with  respect  to  x,  (9)  with  respect 
to  y,  and  substractlng j we  have 

" '^Ixyy  ' ''ixxy  * %yyy  J ‘ 

Which  is  equation  (4. l6)  with  a = 0, 

Thus  the  transport  distribution  for  the  steady  case  is 
precisely  the  same  as  it  is  in  Problem  1»  The  difference  in 
behavior  enters  into  the  non-steady  case  when  the  motion  of  the 
interface  affects  the  motion  of  the  water  in  the  top  layer* 

If  we  set  a = 0,  then  equations  (4,22)  and  (4.23)  are 
the  solutions  for  the  present  Similarly  with  d = 0, 

from  equations  (8)  and  (9)  above 

It  4 ^ 


where  is  given  by  (4,26)*  Then  H2_  may  be  written 


H,  = - 


© + 91/1  + 2ns\ 

©b  ° 

nsX 
b 


(12) 


However,  if  2n^/©b  < 1,  then  may  be  written 

approximately 

e + [9  + H ] 


^2  can  then  be  evaluated  by 


(7) 


If  the  dimensional  constants*  which  were  used  in  Problem 
* The  depth  (D2^-D2)  is 'given "the  same  value  as  D in  Problem  1* 
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f 


1. 


1 

j 


1 are  used  here,  and  if  we  put  b = *005,  then  (12, a)  is  correct 
to  0(10  ),  The  streamlines  and  the  thermocline,  H2  are  shown 

in  Figs,  7 and  8, 

In  Fig.  8 it  can  be  seen  from  the  contour  lines  of  the 
thermocline  that  there  is  not  much  deviation  of  the  thermocline 
from  its  equilibrium  position.  In  particular,  if  the  initial 
depth  be  500  meters,  the  thermocline  does  not  fall  more  than  35 
meters  below  its  average  depth  in  the  southern  half  of  the  ocean. 
In  checking  our  results  with  observation,  we  find  that 
quantitatively  this  result  is  in  poor  agreement  with  field  evi- 
dence, The  definition  of  the  thermocline  in  the  real  ocean  is 
vague,  however,  and  hence  the  two  parameters  Q (corresponding 
to  the  average  thickness  of  the  top  layer)  and  b (the  density 
difference)  are  not  clearly  determined.  In  fact,  they  may  vary 
over  a wide  range  giving  rise  to  a very  considerable  variation 
in  the  deflection  of  the  thermocline# 

In  Fig,  9,  the  vertical  cross  section  of  the  ocean  at 
y^  =0,25  is  shown  for  four  combinations  of  © and  b#  If  we 
consider  the  curve  with  © = 0,0if92  (Dj^  - D2  = 200  m, ) and 
b = 0,0025,  our  result  is  in  good  qualitative  agreement  with 
measurements  of  the  thermocline  off  Chesapeake  Bay  [10],  Quan- 
titatively, the  values  are  out  by  a factor  of  (approximately) 
three# 

Our  solution  shows  a tendency  for  the  thermocline  to 


approach  the  surface  in  the  northern  part  of  the  ocean  (Fig,  8), 
As  a matter  of  fact,  if  © and  b be  chosen  small  enough,  the 
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interface  lies  above  the  free  surface  I Such  a result  is  absurd, 
of  course,  but  the  tendency  of  the  thermocline  to  approach  the 
surface  in  the  northern  part  of  the  ocean  is  clearly  indicated. 
This  fact  agrees  with  observation  since  the  thermocline  actually 
reaches  the  surface  in  the  north, 

Non-Steady  Wind 

In  the  treatment  of  the  non-steady,  two-layer  problem, 
we  shall  neglect  the  terms  with  coefficient  X in  equations  (1), 
(2),  (4),  (5).  For  the  steady  problem,  if  0 and  b are  chosen 
appropriately,  it  has  been  shown  (equation  (12, a)  that  the  error 
Involved  herein  is  small. 

Two  methods  of  attack  have  been  applied  to  the  lineariz- 
ed equations  of  (1)  - (6),  Our  first  procedure  is  that  used 
in  Problem  1,  viz,  a perturbation  in  6 followed  by  a boundary 
layer  analysis. 

The  difficulty  in  the  first  method  of  solution  arises 
from  the  fact  that  the  quantities  with  coefficient  6 are  no 
longer  small,  i. e, , the  magnitude  of  the  terms  is  no  longer 
governed  by  6,  In  particular,  in  the  continuity  equation  (3), 
the  term  on  the  right  hand  side  has  magnitude  6/b  H-j_  (based  on 
the  steady  solution).  In  the  interior  of  the  ocean  where  and 
are  0(1)  and  = 0(0"^),  in  order  for  the  perturbation  in  b 
to  be  valid,  we  must  have  6 < < l/©b.  With  the  dimensional  con- 
stants of  Problem  1,  this  means  6 < < 10”^,  Such  a value  corres- 
ponds to  a wind  period  of  one  hundred  years  or  more. 

If  the  above  results  were  the  only  objection  to  the 
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analysis,  the  problem  as  defined  thus  far  might  still  have  some 
qualitative  value.  Unfortunately,  for  such  a small  value  of  6, 
the  terms  in  the  equations  of  motion  which  involve  a time- 
derivative  become  very  small,  and  we  are  wholly  unjustified  in 
neglecting  the  non-linear  terms  while  still  retaining  these 
time  dependent  terms. 

In  spite  of  these  objections,  the  analysis  for  Problem 
2 by  the  first  method  was  carried  through  but  the  results  were 
not  computed  numerically.  The  analytical  results  are  listed 
in  the  next  few  pages. 


Ui  = ^1  = ^10  ^\l»  % = %0  '’^l 


^2  " ^20  ^^21*  ^2  " ^20  ^''^21’  ^2  “ ^20  ^^21 

where  UgQ  '^20  ~ ^ equation  (49),  U^q,  V^q,  are  given  by 
equations  (3*  22) , (3,  23)and(3,'26)  and  the  remaining  values  are 
given  below, 
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1/3  A 1 (x-r)e"^^^ 

+ e "^A  f e + 


■J 


y sin  nsy  |'[(£  - e^/3)  ^ 
. -1/3 

sin  -)  + [r£ 


r£!^]  ^ =in  . 2£i2 

3^/3  ^ 3 3 


2 ^1/3  „ /t„~1/3  «1/3 


V)3-<=-  (£vl3|_)  , f , . 2.e 

nsy 

■ -1/3 

. vQe  s o _-'-/3 


b©  '3  3 .2  ■-  ‘ 3 

^ sin  i/3e  + r ^sy^^ j,gl/3  ^ 

^ 2 ^“b©"^  3 ^ " 41  cos  • 


• ) e 2 ^ 


^ — y sin  nsy  [(C^  +1/3  C^)  cos  • 


-1/3 


V /o  r~  “1/3  -A 

• ) + (O3  -1/3  C^)sin(.l^ )]  je 


- 2L1 


-1/3' 


V--  = -g  y (y  cos  nsy  + nsy^sln  nsy)(x  - r + e^'^3) 

£ii  9b 

-1/3 


+ a cos  TT 


nsy 

©b 

,2 


sin  nsy  + A2(y)  j.  e 


(x-r) e 


+ a cos  sin  syj-S  cos(.~../3  e j 


I 


/-^“l/3 

2 


o c-1/3  “1/3  1 - M 

+ (2^ 1)  sin  (£il3e )]  ^ 2 

3 v/3  2 ' ^ 


-1/3 


.rr  -1/3 


Xe-i 

-1/3  ^ ^ 


+ a =os  T ]C22Cos(iL!d3|Ji:).C„sln(£is^ 


'32  ' 2 

^21  ~ °©b  ^ “ rx-xe^'^^)  [(2+y^n^s^)cos  nsy+2yns  sin  nsy] 

2/3 


+ a cos  - a cos  sin  nsy  [ (-^~ 

9b  ay  L ©b  J ^ 


-V3 


1/3  /T  1/^  /- 

. iLS^)cos(£i:l|S ) + sln(3L^ 

3 \/3 

+ (el/3  . I)  sln(iy3|±2)+(^  . £)  ^ ^os  • 

^73  2 33 
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r^-1/3  . 2LS. 


-1/3 


e 


- a cos  t: 


9 I sin  nsy  « 


ay 


©b 


-1/3 


* [(\/3^22“  0^2) sin  (—  ^ ) - (^22  ^/3  C -55)003 


. ( 


X \/3e 


-1/3 


- xe 


-1/3 


'32' 


•)]  le 


-1/3 


- a cos  T -i.  - sin  nsy(.^n£*:^ — ^ 

ay  t ©b  3 


+ A^) 


I/3I  (x-r)e 

^ J 

= - ac  cos  y sin  nsyEe^^^^CS-".-^*^-  g 


^ -1/3 

y sin  nsy  [(.^^  - ZLSi- 

Vl 


-2/3  /-  -1/3  2/3 

iL§ S)sin(lZ3e ) ^ (xT£_f 

3 \/3  2 


,-r  -1/3  - 1/3 

osC^LjOi )]e  ^^.±2 


3 \/‘3 

i)  + . 22c^:^  . 

-1/3 

e 1 /-i 


• y sin  nsy  [(Cgg  + \/3 

/,  -1/3  _ -V3  . ■} 

sCiv£^ ) + (C23-- -/3C22)sin(ii-!^^| ) ]e  2 J 

g cos 


T nsy^sin  nsy  1/3 
b©  ® 


^ (x-Ds-1/3 

J e 


(x  - r - r 1 

3 ^ 

+ (1  + ns)  g cos  T nsy3sln 
I 9 >^Q 


1/3 


nsy 
b© 

1/3  /■?  “1/3 

3 \/3  3\/3  2 

1/3  1 /o 


1/3 


(©22  + /3  C 

b©  L 


J ^ -1/3 

~-^)x  ]cos(ii3^3| ; 

/-  -1/3 
I -so)  cos  ^ 


2 

sin  nsy 


-1/3 

- 2LS 

■)  f e 2 


2 
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|[l(x2  + r^)  + (e^/3 


- r)(x  + e^/3)  ] • 


* [ yn^s^  sin  nsy  + ns  cos  nsy  ] - I(cos  + y sin  nsy)[ 

3 ns  J 


g cos  T 


- r)(x  + • 


3 ? 

• [ nsy^sln  nsy  + y cos  nsy  ] + nl  (1  + H§)  • 


• [ ky^  cos  nsy  + nsy3  sin  nsy  - - il«£2£Ji£y]  , 

ns  2 2 

n s 


The  functions  C,  C^,  C^,  C^,  C22J  are  determined 

by  applying  the  boundary  conditions  U]^  = = 0 on 

X = 0,r, 
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1/3  r o o „2 
A sin  nsy  = — i (©n'^s  + 


)S1„  nsy  + LS^JISZ  j 


Cj  sin  nsy  = [(SnSs^  f 2!^)sln  nsy  h-  | 

C(y)  = q(-^  “ -|^(©n^s^  + — ^^^)cos  nsy 

2yns  . ynse^^^  a 

b •'3b©  j 

- “)cos  nsy  j- 

^ j \/3  sin  nsy  J©n^s^  + -i(nsy^ 

2 j © 2 Lb 

- '1  p/^  ^ p p p 


■cos  nsy 


+ Xi- 
© 


+ ~)  \ sin  nsy  + ] (-  4 sin  nsy  + -L. 

ns  J © (^  2 b 2 ^ 

1 5nsy^  ©n^s^\  y *i 

tai  “1^  * ®-g— )sln  nsy  - ^ cos  nsyj 

nsy^sin  nsy  A2  = - e^'^^Cy  cos  nsy  + nsy^  sin  nsy) 
sin  nsy  C22  = - (r  - e^'^^)(y  cos  nsy  + nsy^sin  nsy) 

= - (-^  - [(2  + y^n^s^)cos  nsy  + 2yns  sin  nsy] 

2 2 2 

+ sin  nsy  + cos  nsy  - iL_|-Z_  cos  nsy] 

/-  1/3  p / o 

sin  nsy  Cno  = (-—  - re^'^)  [nsy^sin  nsy  + J— SiSLilSZ] 

2 2 ns 

2/^,y  cos  nsy  , ■^  2 . v 

+ e ^ ^ nsy  sin  nsy) 


2 

+ ^ (nsy  sin  nsy  + 3y  cos 
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The  second  method  of  attack  on  the  non-steady  two-layer 
problem  consists  of  separating  the  expression  for  the  wind- 
stress  into  its  steady  and  periodic  parts,  i.e»,  (1  + a sin  m) 
cos  nsy  = cos  nsy  + a sin  a cos  nsy,  and  treating  each  problem 
separately*  This  method  of  solution  was  also  attempted  in  the 
one-layer  problem.  The  resulting  equations  could  not  be  solved, 
however,  without  recourse  to  numerical  methods.  In  the  present 
case,  we  hope  to  make  use  of  the  smallness  of  the  parameter  b 
in  seeking  a solution. 

In  equation  (3)  the  right  hand  side  may  be  approximated 
by  a/ ax  (-  + aH2)a^  a/3x  (-  dE2/d'V* 

The  steady  problem  with  cos  nsy  as  the  wind-stress  term 
has  been  solved  previously.  For  the  time-dependent  problem,  we 
write 

Hi  = H2=9bH2,  I = Y = . 


Then,  with  the  time -dependent  part  of  the  wind-stress  only,  (1) 
(6)  become 

aa,  al 

ns 6 — - nsyVj^  = - — i + nseAU^  - o sin  cos  nsy 


OVt  -'t 

ns  6 — i + nsylln  = - — i + nseAV, 

ax  -‘-ay  1 


iZi . - 7 

9x  ay  ^ 3x 


aTT  aH  3H 

fS'T  ^ etx  ax 
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8V.  - 3H  3H^ 

b - — + ylJp  = - t + ~-)  + eAVp 

3t  ^ 3y  3y  ^ 


3Up  3Vp  _ 9H 

+ — £ = - Y “ 

3x  3y  ' 3t 


Next  let  us  write  the  wind-stress  as  the  imaginary  part 
of  ae^'''  cos  nsy«  Then  if  we  take  only  the  imaginary  terras  in 
the  remaining  parts  of  the  equation,  the  results  will  be  the 
same  as  those  above. 

Define 


^1,2  = 

The  equations  become 
ins6U'j_  - nsyv^ 
insivj^  + nsyuj^ 

^Ix  '^ly 

16U2  - y^2 

ibvp  + yup 


u 


2x 


+ V 


2y 


u a It 


1,2  “ '^l,2^^>y^  > %,2“ 


ae- 


^h^^2(x,y)< 


3hi 

— * + nseAu- 
3x  1 


3h 


i + nseAV]_ 


= ir  h 


= - -t  ( 


■ 3h]_ 

3h2 

'3x 

,3h;j_ 

+ 

ay' 

3y 

= - iyh 


Y"2  '* 


- cos  nsy 


+ eA^2 

+ eAVg 


The  above  equations  must  be  solved  for  the  six  unknowns. 
The  difficLilty  arises  in  trying  to  match  the  boundary  layer  con- 
tirbution  with  the  Interior  solutioric  To  conserve  space,  we 
shall  not  give  the  entire  analysis  here,  but  shall  confine 
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ourselves  to  the  determination  of  the  boundary  layer  contribu- 
tion and  to  an  indication  of  the  ensuing  difficulties. 

Carry  out  the  following  three  steps: 

l/o 

(a)  Let  X = e i*®**  stretch  x coordinate  near  x = 0. 

(b)  Substitute  2 ~ ^12’  Y ~ 

(c)  Keep  the  leading  terms  of  the  equations. 

The  equations  then  reduce  to 


nsyV]_ 

-1/3 

insie  + nsyU]^ 


= h 


- h3_y  + nsv^^^ 


^ly 


- iyhg 


yvg  = 

4.  -1/3 

15 e V2  + yu2  = 


X (hj^^  + h2^) 

“ ^ ^^ly  ^ ^2y^  ^211, 


^2^  ^2y 


--  - iyh^2  • 


Eliminating  all  the  unknowns  except  h2,  we  find 


where 


^2U^^  - h2^^  - h2^  + L h2  = 0 

b = i5e"^'^^,  L = lyy^Cns  + i)  • 

^ I 


Solutions  are 


ho  = S C.  e 
i=l  ^ 


where  the  Dj^  are  the  roots  of 


4 2 

- D + L = 0, 
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They  are 

%,2 


D 

where 

A 


3,^ 


- (^  5^+2A+2B)  +\[^ ij^+2A+2B)^  2(A+B  - ^ +2) 

2 

^6^+2A+2B  ± fel■^2A+2B)^  - 2(A-fB-  ^ ■ 2) 


- [ i(i+L  + -i)  r+[i  + — L.  - _i]  I . 

L V 2 >-2  + 3 27  r 


The  above  solution  for  li2  must  now  be  substituted  Into 
the  previous  six  equations  and  the  boundary  layer  contributions 
for  u , V h ^ can  be  derived  by  keeping  the  parts  which  ■ 

Xjc- 

0 as  ^ — ^00  • If  the  interior  and  boundary  layer  solutions 
are  added,  the  can  be  evaluated  by  means  of  the  boundary 
conditions  = U2  = V2  = 0 on  x = 0,r, 

Practically,  this  is  an  almost  impossible  task,  and 
numerical  methods  must  be  employed  for  the  whole  procedure# 

In  view  of  this  fact,  nothing  is  gained  by  the  analysis  and  the 
entire  solution  might  as  well  be  carried  out  numerically  from 
the  very  beginning* 
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Since  we  have  been  unable  to  arrive  at  a useful  solu- 
tion for  the  non-steady  ocean  circulation  without  assuming 
negligible  velocities  in  the  bottom  layer,  we  have  no  assurance 
that  our  analysis  is  valid*  Reliable  observational  data  which 
might  guide  us  in  this  matter  are  not  available*  We  may  per- 
haps gain  a little  more  confidence  in  the  results  of  this  inves- 
tigation by  the  following  considerations. 

For  the  formulation  of  Problem  1 it  was  assumed  that 
the  velocities,  and  hence  the  horizontal  pressure  gradient, 
vanish  in  the  bottom  layer*  This,  together  v/ith  the  hydrostatic 
pressure  law,  Immediately  led  to  the  conclusion  that  the  thermo- 
cllne  responds  instantaneously  to  any  motion  of  the  free  surface. 
Naturally,  this  can  hold,  if  at  all,  only  for  sufficiently 
slowly  varying  circulation* 

Some  investigators  are  of  the  opinion  that  the  very 
opposite  situation  actually  exists,  i.e*,  the  thermocllne  re- 
mains essentially  fixed  and  does  not  respond  to  wind  variations 
of,  say,  seasonal  or  annual  periods.  This  is  perhaps  a more 
reasonable  assumption  because  it  is  based  on  the  idea  that  the 
frequency  of  wind  variation  is  much  greater  than  the  important 
frequencies  of  free  oscillations  of  the  bottom  layer. 

Let  us  assume,  therefore,  that  the  shape  of  the  thermo- 
cline  remains  roughly  fixed  in  such  a manner  as  to  result  in  a 
vanishing  time-average  horizontal  pressure  gradient  in  the 
bottom  layer.  That  is  to  say,  the  thermocllne  adjusts  itself 
to  the  mean  wind  distribution  so  as  to  give  zero  pressure 
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gradient  for  the  case  of  a steady  wind  having  this  mean  distri- 
bution# If  we  now  have  a time-dependent  wind,  we  will  have 
non-vanishing  pressure  gradients  in  the  bottom  layer  as  a result 
of  changes  in  the  free  surface  shape.  The  resultant  velocities 
in  the  bottom  layer  will  tend  to  be  uniform  vertically  (except 
as  influenced  by  friction)  provided  the  bottom  layer  has  fairly 
uniform  density  so  that  the  pressure  gradient  is  independent  of 
depth. 

Suppose  we  have  a two-layer  ocean  and  integrate  over  the 
top  layer  only*  If  we  make  use  of  the  assumption  of  a station- 
ary thermocllne,  and  if  the  effect  of  friction  at  the  thermo- 
cline  on  the  transport  in  the  top  layer  is  negligible,  then  the 
resulting  transport  equations  are  essentially  the  same  as  those 
attained  in  Problem  1.  Hence,  the  distribution  of  mass  trans- 
port obtained  in  Problem  1 may  be  expected  to  be  valid  now, 
provided  it  is  interpreted  as  the  distributions  of  transport 
above  the  thermocllne.  Since  this  is  the  transport  usually 
measured,  we  may  still  hope  that  the  results  are  useful, 

6,  Conclusions,  If  the  velocities  in  the  depths  of 
the  ocean  are  negligible,  then  the  horizontal  pressure  gradients 
are  also  negligible  and  the  thermocllne  responds  immediately  to 
a change  in  the  free  surface  height  provided  the  hydrostatic 
pressure  equation  is  valid.  For  such  a case,  the  following 
results  appear  to  be  valid  (within  the  framework  of  subsequent 
approximations  made  in  this  report); 

(i)  For  a varying  wind  with  a period  of  three  months  or 
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more,  the  mass  transport  through  the  Gulf  Stream  responds  to 
the  v/lnd  but  lags  behind  it  at  all  times  except  at  the  instants 
of  extreme  wind  variation  when  the  two  are  in  phase. 

(11)  The  maximum  lag  appears  when  the  wind  is  in  its 
mean  position  and  an  Interval  of  about  nine  days  elapses  between 
the  time  at  which  the  wind  reaches  its  mean  value  and  the  time 
at  which  the  transport  reaches  its  mean  value.  The  actual 
length  of  the  interval,  1, e, , nine  days,  is  independent  of  the 
frequency  of  the  wind  variation. 

(ill)  The  value  of  the  maximum  mass  transport  through  the 
Gulf  Stream  does  not  depend  on  the  frequency  but  only  on  the 
maximum  strength  of  the  wind. 

(Iv)  The  Gulf  Stream  does  not  undergo  any  noticeable  east- 
west  shift  nor  is  its  width  altered  because  of  the  wind  variation. 

For  the  steady  two-layer  problem,  the  streamline  pattern 
coincides  with  that  of  the  one-layer  case.  The  computed  steady 
position  of  the  thermocline  can  be  made  to  agree  qualitatively 
with  the  position  of  the  observed  thermocline  provided  the  two 
parameters  (a)  the  thickness  of  the  top  layer  and  (b)  the  density 
difference,  are  chosen  appropriately. 

At  the  outset  of  our  investigation  we  had  hoped  to  solve 
the  linearl-jed,  non-steady,  two-layer  problem  with  no  a priori 
assumption  concerning  the  vertical  distribution  of  velocity. 
However,  we  were  unsuccessful  in  doing  so  except  for  the  case 
of  a wind  with  a period  of  oscillation  of  100  years  or  more. 

For  such  a low  frequency,  the  retention  of  the  time  derivative 
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terms  in  favor  of  the  non-linear  terms  seems  wholly  unjustified. 
The  only  conclusion  (which  may  not  be  justified  because  of  the 
previous  statement)  resulting  from  this  last  investigation  is 
that  the  transports  in  the  lower  layer  are  of  the  same  order  of 
magnitude  as  the  out-of-phase  transports  of  the  upper  layer. 

In  view  of  the  statements  made  at  the  end  of  Section 
the  results  listed  for  the  one-layer  problem  are  approximately 
valid  for  the  non-steady  two-layer  problem  provided: 

(a)  The  therraacllne  adjusts  itself  to  the  mean  wind  dis- 
tribution and  remains  fixed, 

(b'  The  mass  transports  of  Problem  1 are  Interpreted  as 
the  transports  in  the  upper  layer. 

The  assumption  of  hydrostatic  pressure  is  not  necessary 
for  the  solution  of  the  mass  transports  in  the  steady  problem. 
Wherever  the  results  of  this  analysis  permit  a compari- 
son with  observation,  good  qualitative  agreement  is  achieved, 
but  the  quantitative  results  are  off  by  a factor  of  about  three. 
In  view  of  the  many  idealizing  assumptions  made,  however,  no  more 
than  qualitative  agreement  could  be  hoped  for, 

A number  of  features  have  been  left  out  of  the  present 
model.  Changing  topography,  non-linear  terms,  variable  eddy 
viscosity  and  many  other  featvires  could  combine  to  change  the 
results  noticeably.  However,  the  analysis  of  the  problem  in- 
cluding most  of  the  features  which  were  omitted  in  our  model 
would  probably  require  a numerical  treatment. 
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Appendix  1 . Transformation  of  the  Differential  Equations  from 
Spherical  to  Rectangular  Coordinates. 

Consider  a rotating  spherical  coordinate  system;  let 
r be  the  radial  distance  from  the  center  of  the  sphere,  0 the 
colatitude,  «p  the  meridlanal  angle.  The  equations  of  motion  are* 


^ - Q^r  sin^©  - 2wQ  sin  © = -lA^-g'  + l(v»AiV) 


p 8r 


u 


+ u ^ + x ^ + _JL_ 9v  ^ uv  _ y.ic.qt_Q  + gin  Q nog  Q 

9t  0r  r 9©  r sin  © a<p  r r 


- 2wQ  cos  © = - i i + i(V»  Ai  V)  V 

p r 9©  p 


^ j.  ,,  9w  j_  V M 


+ + ¥ — s + 2r  Q cos  © 

0t  9r  r 0©  r sin  © 09  r r 


+ 2uQsln  © = - 1 1 — _ + i(V  • A4  v)w 

p r sin  0 9cp  p ^ 1 

where  ^ is  the  material  derivative  of  the  radial  velocity  in 

JJX 

terms  of  spherical  coordinates 


g'  denotes  the  gravitational  force 

- V * Ai  V = AV  ^ + T a—  (Ao  -^)  and  denotes  the 
p p0rJ9r 

Laplacian  operator  for  the  two  dimensions  © and  9 , 

V/e  shall  neglect  the  radial  acceleration  and  shear  terms 

arising  as  a result  of  the  velocities  relative  to  the  rotating 


We  shall  not  consider  the  non-linear  terms  or  the  viscous 
terms  in  the  radial  equation  of  motion;  hence,  this  equation 
is  written  in  operator  form  only. 
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sphere.  We  then  have 


1 ^ = 
p 3r 


g 


(1) 


^ + V 3v  -I.  uv 

3t  9r  r 90  r sin  © 9<p  r 


2 

w cot 


0 


sin  0 cos  6 


-2w  Q cos  0 = - ii^  + Ay^v+i  ™-(Ao 
p r 9©  p 9r  3 


mm 

9r 


9v  ^ + V w 9w  ^ vu 

9t  9r  r 9Q  r sin  © 9<p  r 


yW-^qt,Q  + 2vQcos  © 
r 


- 2u  Q sin  © = -i 


^4 — s ^ + Ay^w  + i A (A,  -1^) 

p r sin  9 9cp  p 9r  3 9r 


(3) 


vhere  g = g’  - apparent  gravitational 

force.  The  viscous  terms  for  equations  (2)  and  (3)  are 


A y^v  = A. 


cot  © ^ + ^ + -X  - . 2_^os^  ^ 

Q®  9©^  sin^©  9<p2  sin^©  sin^© 


Ay^w  = A- 


{ 


cot  © 


9v 

9© 


9^w 

9©^ 


sin*^© 


w 


w 


09 


sin^© 


2 cos  © QV 
sin^©  99 


Since  the  region  of  interest  to  us  consists  of  a very 
thin  layer  on  the  surface  of  the  globe,  we  shall  approximate 
r by  R,  the  mean  radius  of  the  earth,  whenever  r appears  in 
undifferentiated  form.  At  the  same  time  let  us  define  a new 
east-west  coordinate  by  x = yR  sin  ©,  a north-south  coordinate 
by  y = R(2  - ©)  and  a vertical  coordinate  by  z = r.  Then 
equations  (l)-(3)  become 


1 ^ _ 
p 9z  ~ 


g 


(^) 
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3v  j.  „3v  „ 3v  , 3v  . uv  w^cot  © . o2_  _4„  ^ 

R — ^ ® ® 

r 

-2wQcos©=14E  + aJ  - coO  ^ + a£v  ^ ^ 

P ^ [_  R 3y  ^ ^ 


- _ V ^ 

R^sin^Q 

- ^ Qot  Q ^ "V+  1 A.  CAo  (^) 

R 3x  J p 3z  ^^3  Qz^ 

9v  . ,-^  . ,,8w  . wu  uw  cot  © . o Q o„  n«^v^  o 

3t  '*'  " ^3y  ^ ^'dx  R “ R + 2v  w cos  © - 2u  S2sin  © 


= - i + A 

p 3x 


cot  © ^ 

3y 


R 


<j  w I 8^w  1 3^w 

t + ' • -A'  T 

3y2 


3x 


w 

R^sin^© 


2 cot  © 3v 
R “■"  d~ 


3jr  L 

J 


+ 1 ±.  (Ao 

p 3z  ^ 3 3:r 


(6) 


Since  R is  very  large,  we  shall  neglect  terms  divided 
by  R.  We  can  do  this  provided  the  region  is  sufficiently  far 
removed  from  the  poles  (©  = 0,x)  where  cot  © becomes  infinite. 
The  velocity  component  u is  assumed  to  be  much  smaller  than  the 
components  v and  w so  that  we  can  neglect  u throughout  the 
equations  of  motion. 

Ordinarily,  one  uses  the  velocity  components  u,v,w  to 
correspond  to  the  directions  x,y,z  respectively.  In  equations 
(4) -(6),  u,-v,w  correspond  to  z,y,x  respectively.  The  negative 
sign  was  carried  over  from  the  definition  of  v which  was  defined 
positive  southward.  If  we  revert  to  the  more  familiar  nota- 
tions and  vjrite  u'  = w,  v'  = -v,  w'  = u,  we  have  for  equations 

(4) -(6) (with  the  terms  with  coefficient  1 and  all  terms  con- 

R 

taining  w’  neglected) 
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|u'+  u'^’+  v'|B'  - 2Qv*  sin  a) 
at  ax  3y  R 


= . 1 ^ + A |^'+  2!^ 

pax  I 9y2 


+ 1 ^ 
p az 


(7) 


V'#'  + 2Qu»  sin  sin  © cos  © 

at  ax  ay  r 


= . i 2£  + ii' 

p ay  ax^  ay'^ 


D ai  ^ 3 aT^ 


(8) 


1 ^ = - 

p 8z 


g 


(9) 


If  the  above  procedure  be  carried  out  for  the  con- 
tinuity equation,  the  latter  becomes 


au»  ^ = 0 (10) 

ax  ay  3z 


In  making  the  transformation  from  spherical  to 
rectangular  coordinates,  we  must  consider  the  distortion  of 
the  spherical  surface  as  a result  of  the  mapping  process. 
Specifically,  a rectangle  in  the  rectangular  system  maps  into 
a region  on  the  sphere  in  such  a manner  that  the  east-west 
distance  remains  constant  and  the  right  angles  between  the 
lines  X = const,  and  y = const,  map  into  obtuse  angles  be- 
tween the  lines  on  the  sphere  corresponding  to  x = const,  and 
y = const.  Thus,  the  mapping  is  not  conformal. 

With  the  above  transformation  we  have  mapped  a 
spherical  surface  onto  the  plane.  Our  real  aim,  however,  is 
to  map  the  equilibrium  surface  which  establishes  itself  as  a 
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result  of  the  interaction  of  centripetal  acceleration  and 
gravity,  onto  the  plane.  We  shall,  therefore,  neglect  the 
small  difference  between  the  true  equilibrium  surface  and 
the  sphere. 

The  apparent  gravity,  g,  in  (1)  acts  perpendicular 
to  the  spherical  surface.  V/e  shall  now  consider  g to  act 
perpendicular  to  the  equilibrium  surface.  We  must  then 
drop  the  term  Q^r  sin  © cos  9 from  the  6 equation  since,  in 
reality,  this  force  combines  with  g acting  normal  to  the 
spherical  surface,  to  give  rise  to  a resultant  normal  to  the 
equilibrium  surface.  Finally,  g = g'  - Q^r^sin^©)  is 

assumed  constant.  The  final  result  of  the  approxihaate  trans- 
formation is  to  map  the  equilibrium  free  surface  of  the 
ocean  onto  the  x-y  plane,  with  the  apparent  force  of  gravity 
acting  normal  to  this  plane. 


All-101 


8^. 


Appendix  2.  T^epject  of  the  Hon-Llneer  Terms. 


Jonslder  the  Integrated  equations  of  inotion  of  section 


3. 


-h 


3u' 

■al- 


dz  + 


u 


-h 


' dz  + V*  ~i*dz  - py  v'dz 

ax  \ 3y  , 

u -h  (J  -n 


lau' 

3t 

-h 


dz  + 


= _ g D + A 

ax 


Au'  dz  + Tx 


-h 


u 


• 9v'  dz 
3x 


-h 


+ v'  "l^'dz  + P,y  r ^ u*dz 

J -h  J -h 


sd|^  + a|'’^Av*  dz+'^y 


(1) 


(2) 


where  v'e  have  linearized  the  pressure  term  in  accordance  vrith 

remarks  to  be  nade  later  in  sections  3 £>•'1^'-  5»  are  now 

X y 

the  wind-stress  components  of  section  3 divided  by  p.  Assume 


,kz  _ r.r.,  „ 4.\  J<z 


u'  = u(x,y,t)  e , V*  = v(x,y,t)  e‘' 


i.G,  the  velocities  decay  exponent ia  11, y v/ith  depth. 


Then, 


au  .kz 
at  ® 


-h 


+ i u 

2 ax 


2 3y 

•h 


T) 

-h 


2Q  sin(^)  V e 
R 


kz 


-h 


= -gD|Jlk.  + AAu  e’^^ 
3x 


n 

-h 


+ Tj.k  ^3) 
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w_  y 


^ ] 

+ ^ 

u 2V 

'>  + 1 V p 

C 

-h 

; ax 

-h  ^ 

-h 


+ 2Q  sin(^)  u 
R 


n = - 

-h 


k + Mv 


+ T k (4) 

y 

-h 


Approxi:;iato  the  exponentials  at  their  limits  by 
« 1,  e“"^  23  0,  Q?hen  (3)  and  (4)  become 

V = - S D |ak  + Mu  + (5) 

« = - g D |Sk  + MV  + vyk  (6) 

Linearize  the  Coriolis  pareraeter  by  2Q  sin(^)  ~ 

j.\  «y 

where  p = ^ , Taking  the  derivative  of  (6)  with  respect  to 
X and  the  derivative  of  (5)  with  respect  to  y and  subtracting, 
we  have 


_2.  (2v  _ 2Uw  1 

at  ax  ay  ^ 


^ + u + §2  ^ + ~9^v  _ ^ 9u  ^ ^a^u 

3x  ax  5x2  ax  ay  axay  “ ax  ay  axay 


-§ 


- v-^ 


y oy 


^y^”5x 


ay'  '3x  ay' 


+ k(^  - ^ 
^ax  ay 


) 


(7) 


Choose  V = 0 , T = - (W’  + r'  sin  wt)  cos  ny. 
y X 

We  shall  non-dimensionalize  the  velocities  so  that  they 
are  of  order  unity  in  the  interior  of  the  ocean.  It  is  con- 
venient* to  choose 


% 

The  choice  of  the  non-dimensional  quantities  is  motivated  in 
section  4. 
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sx' , y = sy' , T 


W r * 


Y = 


kW 


? 2 

H = 

V/ 


w = nW  , r 


n r» 


In  this  notation  and  with  the  prescribed  form  for  x ^ and  Ty  , 
equation  (7)  becomes 


+ 


Y 


^ M 

.ax'  ax' 


M .21 

0X' ay' 


+ V- 


2 

a V 


ax'3y' 


. au  M - U M su  - v-ain 

ax'ay'  ax'ay'  ay'ay'  gyi2 


+ y 


au  + M 

ax'  ay 


+ V =eA' 


21  - 
ax'  ay' 


- (1  + a sin  T)sin  nsy' 

(8) 


The  integrated,  non-dimensionalized  continuity  equation  be- 
comes 


au  + a^  = - ^ M 

ax'  ay'  ax 


(9) 


If  we  expand  the  velocities  and  the  height,  H,  in  a 
series  in  6 , then  the  solution  can  be  looked  upon  as  the  sum 
of  a quasi-steady  part  plus  a number  of  out-of-phase  contribu- 
tions. If  b is  small  enough  we  may  be  justified  in  keeping 
only  the  first  two  terms  of  such  a series  as  a fairly  accurate 
representation  of  the  complete  series. 

Hence,  let 

U = Uq  + 6 U]_  + 6 ^U2  + ...  , V = Vq  + 6 Vi  + 6 ^V2  + ... 

H = Hq  + 6 + 6 ^H2  + ... 


u 


— A 
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Y 


au.  av. 


ax'  ax' 


- Vr 


a^Uo,  av^ 

On  ,Tr  -A  r'^O 


9y'' 


■]  + Vo  = eA  [: 


au, 


? ] 


■8x'  8y' 

- (1  + a sin  T)sin  nsy' 


au^  av^ 

— 2 + — 2 = 0 
ax'  ay» 


(10) 


(11) 


The  first  order  equations  in  6 are: 

a /av^  au_^  . ail  av^  au,  av^  a^v, 

aT  ax'*  " + Y [ sr:7r  aTt  ~-7r  p + 


2v  q2v 


ay' 


ax'  ax' 


ax'  ax' 

9H, 


ax 


,2 


ax 


i2 


+ ...]  - y — 2 + V]^  =e  A'  (—7 7)  (12) 

ax  ax' 


ay 


a^  ^ av^  _ ^ 


(13) 


ax'  ay'  aT 

Munk,  Groves,  and  Carrier  [?]  have  shown  that  the  effect 
of • the  non-linear  terms  in  [10  ] is  quantitative  and  that  these 
non-linear  terms  can  be  neglected  as  compared  to  the  Coriolis 
term,  Vq.  The  relationship  of  the  non-linear  terms  to  the 
Coriolis  term  in  equation  (12)  is  essentially  the  same  as  that 
in  equation  (10),  This  fact  can  be  shown  by  considerations 
based  on  orders  of  magnitude.  Wo  choose  a typical  non-linear 

auo  aVo  , , au.  av, 

term  in  each  equation,  y Y In 

(12),  and  compare  it  to  the  Coriolis  terms  in  that  equation, 

Vo  in  (10)  and  Vj^  in  (12). 

In  the  solution  it  is  shown  that  Uq,  Vq,  U]_,  V]_  and  all 
their  derivatives  are  of  order  unity  in  the  interior  of  the 


Equation  (10)  with  a = 0 is  the  same  as  that  of  Munk,  Groves, 
Carrier  [7 J . 


All-101 


88 


ocean.  Near  the  boundary  x*  = 0,  it  is  shown  that  Uq  = 0(1), 
Vq  = 0(e"^'^^),  = 0(e“^^3),  Y-^  ~ 0(e~^/3)  and  has  the 

effect  of  multiplying  the  magnitude  of  a term  by  0(e“l/3), 
Based  on  these  results  the  terms  to  be  compared  are 
given  in  the  table  below. 


Interior 
Vq  = 0(1) 


Y ^ = yO(l) 
' ax'  ax'  ' 


= 0(1) 


Near  x'  = 0 
Vo  = 0(e"l/3) 


8U^  av_ 

^ aSc'’  a3c^  ~ ^ 

= 0(e“2/3) 

8Uo  8V-,  ^ 

Y — 2 — £ = Y0(e  ^'•^) 

' ax'  ax'  ' ^ 


Thus,  in  the  interior  in  each  case  wo  have 0(1)  vs, 
Y0(1)  , Near  the  boundary  x'  = 0,  in  each  case  we  must  com- 
pare 0(1)  vs.  Y 0(e~^'^3) , Hence,  the  relationship  of  the  non- 
linear terms  to  the  Coriolis  terms  is  essentially  the  same 
in  the  tv/o  sots  of  equations.  It  would  seem  therefore  that, 
if  the  non-linear  terms  can  bo  noglocted  in  the  steady  equation 
(10),  they  can  also  bo  neglected  in  the  first-order,  non-steady 
equation  (12), 
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Appendix  3 . Hydrostatic  Pressure  Assumption . 


The  results  In  the  main  body  of  the  report  are  based  on 
the  assumption  that  the  vertical  equation  of  motion  can  be 
approximated  by  the  hydrostatic  pressure  equation.  Although 
this  approximation  Is  probably  sufficiently  accurate  for  the 
problem  under  consideration,  It  may  warrant  a few  further 
remarks. 

Consider  the  steady,  linearized  problem.  The  equations 
of  motion  with  a linearized  Coriolis  term  are 

-pyv.  =-l||  + AAu’  (1) 

pyu.  = - i + MV  + i ^(Aj  II'  ) (2) 


and  the  continuity  equation  is 

+ M&v LI  ^ 3.,Cpw'.).  _ 0 

3x  ^ 3y  ^ 3z  ^ 


(3) 


Equations  (1)  and  (2)  can  be  multiplied  by  the  density 
to  yield 


- = - li  ^ A*  ^ 


where  we  have  written  AA(pu')  for  ApAa'  and  AA  (pv')  for 
ApAV . This  approximation  is  certainly  permissible  since 
these  terms  represent,  in  the  first  instance,  only  very  rough 
approximations  to  the  true  state  of  affairs  in  turbulent  motion. 
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If  we  integrate  (l)-(3)  from  a depth  z = - h(x,y,t) 
where  the  motion  is  assumed  negligible  to  the  free  surface 
z = n (x,y,t) , then 


dz  + AAU  + 

\ ax 

J-h 


+ Py  U = - 


^ ^ dz  + AAV  + A^  1^’ 
U-h 


-h 

-h 


0 


(^) 

(5) 

(6) 


where  the  non-llnoar  terms  resulting  from  the  interchange  of 
derivatives  and  integrals  in  the  viscous  terms  have  been 
neglected. 


The  terms 


j az 


n 

-h 


= Tjc 


i 3v' 

3 ai 


= T., 


-h 


provide  the  wind-stress  components  at  the  free  surface  (sob 
sec.  3 of  report).  The  depth  z = -h  has  been  chosen  as  that 
depth  where  the  velocities  are  negligible  so  that  the  contribu 
tions  of  the  above  terms  at  the  lower  limit  are  negligible. 
When  the  term  in  the  continuity  equation  is  integrated, 
it  provides  a contribution  involving  a time-derivative,  viz., 

, so  that  it  vanishes  in  the  present  problem. 

The  pressure  terms  are 


f 1 8e  dz  = _3. 

8y  ay 

J-h 


n 

p az  - |-n 

- 

3x 

P-h 

-h 

n 

p dz  - ^ p 

- ^ 

p , 

ay  n 

-h 
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where  is  p evaluated  at  z = t]  , and  p ^ is  p evaluated  at 
z = -h. 


If  the  free  surface  be  considered  a surface  of  zero 
pressure,  then  p^  = 0. 

Defining 


P 


An 

p dz 

J-h 


we  have  for  equation  (4)  and  (5) 


^ - i " i p-h  * " s (7) 

^ - f " i p-h " ^ V 


A stream  function  ij)  can  be  defined  byU  = -.^,V  = +|^ 

so  that  (6)  is  satisfied  identically.  Taking  the  derivative  of 
(7)  with  respect  to  y and  (8)  with  respect  to  x and  subtracting, 
we  obtain 


A AA  ^1)  - = II 


9P-h 

3y 


-&T 


-Jj 

ax 


(9) 


Since  z = -h  is  the  depth  where  the  velocities  are 
negligible,  the  third  equation  of  motion  below  this  depth  re- 
duces to  the  hydrostatic  pressure  equation,  - = gp  > if  p 

0 z 


is  constant  along  z = -h(x,y,t>.  Then  = gp~  > 


cy 


9y 


9P-h 

3x 


= 6P-i^  • Vfith  these  results  substituted  into  (9),  we  have 
ox 


(10) 


All-101 


92 


! 


If  boundary  conditions  are  imposed  and  if  t and  t 

y ^ 

are  specified,  the  problem  defined  by  (10)  can  be  solved  (see 
Appendix  5).  Thus  for  the  analysis  of  the  steady  state  prob- 
lem, the  only  necessary  assumption  concerning  the  pressure  and 
the  density  is  that  the  density  be  constant  along  the  surface 
below  which  the  velocities  are  negligible. 

If  the  height  z = -h  is  approximated  by  a constant, 

«■ 

then  the  derivatives  of  the  pressure  terms  in  (9)  vanish  and 
no  assumption  need  be  made  concerning  the  density  along  the 
surface  z = -h. 
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. 4(a) . Derivation  of  Re 


sen  T and  ri  . 


With  the  density  distribution  given  by 


P = P 


T)  > z > T 


p = p Q [ 1 + c(T-z)  ] T>z>T-d 


P=P_h=PoI^l  + G<i^  T-d>z 


we  can  find  a relationship  between  T and  by  considering  the 


conditions 


1 ^ 

p ax 


= 0 , i t = 0. 

z=-h  P z=-h 


The  hydrostatic  pressure  equation  is 


1 M 

p OX 


p = g pdC 


P z=-h  ^ l)-h  ^ 


= 0 T]  2 z > T 


p^c  T > z > T-d 

o ax  — 


’-d  > z 


Hence, 
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9^ 


i ^ 

p 8x 


z=-h  P-h 


Po° 


T-d 


|T  dC  + £o  M 
3x  p 9x 


= 0 


or 


od  |T  + |a  = 0 

9x  9x 


9T  _ - ^ 

ax  “ cd  ax  --Kp  9x 


Where  Ap  = p_^  - p^  . 
Similarly, 


aT  _ L2  M 
"Sy  “ “ Ap  "Sy 

Integrating  (1)  and  (2),  we  have  finally 


Pq  di 


T = - 


Po 

Sp 


ri  - C 


where  T = - C when  t]  = 0. 


(1) 


(2) 


Appendix  4(h) . Derivation  of  Integrated  Pressure  Terms , . 

In  order  to  compute  the  terms  f ^ i ^ dz,  P ^ 1 ^ dz, 

(J-h  P 9x  ’ J -h  p 9y  ’ 

we  must  divide  the  region  of  integration  into  three  separate 
parts,  viz., 

n T)  pT-d  pT  p Ti 

J -h  j -h  J T-d  J T 


P T OT-d  ftTl  nT  nT) 

dz  = g 1[  |AdC]dz  + g i[  |fidC]dz 

J -h^  9^  J -h  P J z J T -d  P y z 9x 


(1) 
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Using  the  values  of  ^ for  the  three  layers  listed  in  Appendix 
4(a),  we  have 

= p^c  (T-z)  T-d  < z < T 
= 0 T < z 


T hen 


cd+1 


gp  ^ j ^ i dz  = g T-^-T  [T-d+h]  - g ^ i log(--j-  ■■) 
ax  (J_vi  p ^ 1+cd  ax  c cd+1 


+ g |5  (t]-T). 
ox 


Let  us  put  those  values  into  (1)  and  at  the  same  time 


use 


8T  P-o  ^ _ 1 dp 

ax  ^ - ST  8x  ~ cd  8x 


f '-^1-  c 


''  i ^ dz  = 

-h  P 


g T|-i-T  [T-d+h]  + g rp—  [T-d+h] 

8x  1+cd  8x  1+cd 


- £ Qa  - g ^ 1 |Q  log(  1 ) - g 1 log(-J~-) 

c 8x  S^p  c 8x  cd+1  c 8x  cd+1 
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+ + g(^  r,  + D)  M 

ax  ^ ET  ' dx  ^ ^ dx  Ap  C ^ 


But  , 


log  ^ = - log  ^ = - log  = 


Since  the  term  ^ is  small  wo  can  v/rito 
•■^o 


logd  + ^£) 

r A 


^ ^ (.46)2 

O *^0 


Henco  (2)  becomos 


f_n  ^ dz  = g (c . ^ 1^ 


Similarly, 


i ^ dz 

-h  p 9y 


= g (C  + ^ 


W ^ T]  + 4) 

W ' ? 3y 

= g (D  + ^ T)  ) ^ 
Ap  ' ay 


where  D = C + ^ 
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Appendix  . An  /Example  of  Boundary  Layer  Technique . ♦ 

In  this  section  we  shall  discuss  the  application  of 
the  boundary  layer  technique  to  the  solution  of  the  problem 
defined  by  the  equation 

eAAi})  - = (1+a  sin  t ) sin  nsy  (1) 

and  the  boundary  conditions 

\|)  = = 0 on  X = 0,r  (2) 

^ = 'I’yy  = 0 y = 

The  nature  of  the  boundary  layer  problem  Is  characterized 
by  throe  features;  (1)  the  problem  Is  non-dlmenslonallzod  so 
that  the  size  of  the  domain  has  lengths  of  order  unity | (2) 

the  coefficient  of  the  most  highly  differentiated  term  is 
small  compared  to  unity;  (3)  the  remaining  terms  have  coeffi- 
ci'^nts  of  order  unity.  The  problem  to  be  considered  here  has 
already  been  put  into  a suitable  non-dimensional  form. 

If  i])  were  everywhere  a smooth**  function  of  its  arguments 
and  of  order  unity,  then  It  should  be  possible  to  determine  a 
good  approximation  to  the  solution  by  neglecting  the  term  with 
coefficient  e(e<<l)  and  by  considering  the  remaining  equation 


* For  an  interesting  account  of  boundary  layer  technique,  in- 
cluding the  treatment  of  non-linear  problems,  the  reader  is 
referred  to  [8]  . 

By"smooth"  wo  mean  that  ;|)  has  no  large  derivatives,  i.e. , 

» 'l^x’  'J'xxxx’  order  of  magnitude. 
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“ (1+0^  sin  t:.)  sin  nsy  (3) 

Thus,  a possible  solution  is 

= (1+a  slnx)  sin  nsy  [ - x + Cix(y,T)  ].  (4) 

We  arc  now  faced  with  a dllemna,  however.  as  given 
in  (4)  provides  one  arbitrary  function  of  y and  t to  satisfy 
the  four  conditions  on  the  boundaries  x = 0,  x = r.  If  our 
assumption  that  \)j  is  evorywhoro  a smooth  function  is  correct, 
then  we  are  at  a loss  to  find  a complete  answer  to  the  problem. 
For  if  ij)  and  its  derivatives  have  the  same  order  of  magnitude 
everywhere,  the  only  possible  solution  is  of  the  form 

+ 0(e)  and  it  is  not  possible  to  satisfy  all  boundary  con- 
ditions. 

It  is  obvious,  therefore,  that  \|)  cannot  be  smooth 
everywhere.  In  particular,  in  order  for  the  full  solution 
to  be  different  from^|jj_  + 0(e),  at  least  one  of  the  terms, 
^xxxx’  '^xxyy’  ^yyyy  order  e“l  in  some  part  of 

the  domain  under  consideration  so  that  the  approximation  of 
neglecting  terms  of  order  e will  not  reduce  the  order  of  the 
differential  equation.  If  \1>  is  smooth  away  from  the  boundaries 
and  if  derivatives  with  respect  to  x are  large,  so  that 

is  of  order  e"l,  near  x = 0,r,  then  the  problem  is  one 
of  the  boundary  layer  type.  We  shall  proceed  formally  on  the 
assumption  that  this  is  true,  realizing  that  if  it  is  not  the 
case,  we  shall  be  led  to  a contradiction. 
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The  solution  may  now  be  written  as  the  sum  of  two  parts- 
given  by  (4)  (the  "interior  solution"),  being  sensibly 
large  only  near  the  boundary  and  negligibly  small  in  the  in- 
terior, (the  "boundary  layer  contribution").  VJe  must  now  try 
to  determine  the  boundary  layer  contribution. 

The  nature  of  the  total  solution  Itself  is  the  guiding 
factor  in  the  investigation.  We  have  supposed  that  near  the 
boundaries  x = 0,r,  \()|^  has  largo  derivatives  with  respect  to  x 
while  '4ij_  is  over5rwhere  smooth  and  of  order  unity.  Thus,  if  we 
write  our  solution  in  two  parts,  i.e.,  "the  differential 

equation  can  be  written  in  the  form 

eAAt|)j_  + eAA  = (l+a  sin  m)sin  nsy. 

Now  the  term  eAAil)^^  is  of  order  e,  the  terms  underlined  twice 
are  of  order  unity  and  the  order  of  magnitude  of  the  terms 
underlined  once  is  as  yet  undetermined.  Since  the  terms  in 

are  to  have  derivatives  with  respect  to  x which  are  (assumed) 
large,  we  have  » 1.  Hence  at  least  one  of  the  terms  of 
eAA^J)|^  must  bo  as  large  as  in  order  to  balance  this  term. 

The  equation  will  then  be  satisfied  approximately  if  wo  write 


and 


- = (l+a  sin  t)  sin  nsy 


"“'I’b  - ♦bx  = ° 


C5) 


We  must  now  integrate  these  equations  and  then  add  the  two 
solutions  and  to  form  the  comploto  solution  \()  , 
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The  solution  to  the  first  of  the  two  equations  is  given 
by  (4).  Since  the  complete  solution  will  only  be  approximate, 
in  that  terms  of  order  e have  already  been  neglected,  ^i)^  need 
only  be  determined  approximately. 

It  is  suggested  by  the  above  considerations  that  we 
find  a formal  method  for  writing  our  equation  so  that  the 
magnitudes  of  the  terms  are  expressed  by  the  coefficients  and 
that  the  derivatives,  etc.,  be  of  order  unity.  VJc  can  do  this 
by  stretching  the  x coordinate  near  the  boundary  i.e.,  by  de- 
fining a now  X coordinate  so  that  a particular  distance  in  x 
becomes  a much  larger  distance  in  the  now  coordinate. 

Formally,  wo  operate  as  follows.  Let  x bo  replaced  by 
the  coordinate  ^ such  that 

X = 

where  n is  to  be  determined.  Then  the  equation  (5)  becomes 

yy  + e^,t>yyyy  = 0. 

In  choosing  n we  note  that  it  must  be  positive  if  the 
X coordinate  is  to  bo  stretched.  Thus  of  the  terms  which 
originally  had  coefficient  e,  is  the  largest 

since  it  has  the  largest  coefficient  (n.b. 

^b^^yy’  '^byyyy  order  of  magnitude).  This  term 

is  matched  with  ^b^’  remaining  large  term  in  the 
differential  equation,  and  by  equating  the  coefficients  of  the 
above  two  terms,  we  have  n = 1/3. 
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Thus  we  get 

^ ’•’bSCyy  + ^ Vyyy  - ' ° 

or 

■ '^hK  " 

Now  if  be  expanded  into  an  asymptotic  series  and  if 
we  keep  only  the  first  term  in  the  series  (for  all  practical 
purposes,  this  amounts  to  neglecting  the  0(e~/3)  terms),  we 
have 


■J)  - il)  =0 


(6) 


The  solution  to  (6)  is 


2xi 


'•'h  = + C22(y>'>^)e^  + C^p(y,T)e^®  3 


32 


,ki 


We  have  specified  that  this  solution  is  to  become 
negligibly  small  as  the  distance  from  the  boundary  increases. 
Thus  letting  ^-ioo,  we  note  that  it  is  necessary  that  0-^2  = 
C22  = 0 since  neither  0^2  tends  to  zero.  Hence,  for 

the  region  near  x = 0,  we  have 

r p 2x1  r p 4xi 

'I'b  = 3 + 

1/3 

or,  changing  our  coordinates  back  to  x by  means  of  x = e ^ , 


All-101 


102 


xe-1/3  e 


3 + C.  „(y,T:)e 


xe-1/3  e 


^icl 


'b  " ''32'-'’“'''  ' "^2 

For  the  boundary  near  x = r,  we  now  define  ^ by 


(x-r)  = e’^^ 


and  specify  that  the  solution  vanish  as  ^-oo,  i.e.,  as  the 
distance  into  the  interior  part  of  the  ocean  increases.  By 
a similar  analysis,  we  find  that  near  x = r, 

r 

% = Ci3^y»'^)  + C^^(y,'c)e^  + C^^(y,T;)e^®  3 

le  1+Si 

+ 3 

In  order  for  to  tend  to  zero  as  - oo  , it  is 
necessary  that  C]_3  = C33  = = 0.  Hence 

% = G23(y,T^)e^  = C2^(y,x)e^ 


Hhe  total  boundary  layer  solution  can  be  written 
(x-r) 


= CgCy,  T)e 


+ G^Cy,  T)e 


xe-1/3 


+ Gi^(y,  T;)e 


xe 


-1/3 


(7) 


The  solution  throughout  the  domain  consists  of  (4)  and 


(7)  or 


rjj  = = (1+  a sin  t )sin  nsy|^[  - x + C^(y,  t)  ] 

»-V3e>td  I 
+ C^.(y,T)e^  3 J 


(8) 
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An  application  of  the  boundary  conditions,  \j)  = ■ij)  = 0 
on  X = 0,r,  yields 


= (1+  a sin  T ) sin  nsy 


-x+r-e 


1/3 


1/3  (x-r)e~^/3 

+ e e 

" ' 2 


+ [(e’-^^.rjoos  t ((/3  eV3  . ^)sin  (x}!!^^] 

1 ^ V3 ^ , 

3 ' 1 

^ I } 

The  term  1 is  valid  throughout  the  ocean.  Near  x = 0,  3 t>e- 
comes  as  important  as  1 and  gets  negligibly  small  as  x in- 
creases. Near  x = r,  2 and  1 together  form  the  solution  but 
2 tends  to  zero  as  x decreases. 

Perhaps  a few  remarks  should  be  made  as  to  the  specific 
choice  of  sin  nsy  for  the  total  y dependence  of  the  solution. 
The  particular  choice  of  sin  nsy  satisfies  the  boundary  con- 
ditions 'l>=i)  =0ony  = 0,  y = l,  and  is  supported  by  the 

V V 

specified  wind  distribution.  Thus  we  were  not  forced  to  resort 
to  a boundary  layer  analysis  to  satisfy  the  four  boundary 
conditions.  Of  course,  such  a simple  choice  is  not  always 
possible,  and  one  might  have  to  resort  to  methods  for  refining 
the  interior  solution  in  other  problems  in  order  to  satisfy 
the  necessary  boundary  conditions. 


The  rectangular  oeeah  with  dirhensions  in  both  the 
dimensional  and  the  non-dimensional  coordinate  system. 


Fig.  3.  i\lon-dimensional  mass  transport  function  V near  x'-  r (eastern  boundary). 


one-layer  and  two-loyer  oce 


Fig.  8.  Contour  lines  of  the'  thermccline  for  the  steady  problerri.  Numbers  denote  meters,  above  equiiibium  position. 


